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P  RE7:  ACE 


The  computer  program  (CURVES)  described  in  this  Memorandum  was 
developed  in  support  of  estimating-relationship  research  e.  orts  being 
conducted  in  the  RAND  Cost  Analysis  Department.  This  program  represents 
a  compilation  of  various  parts  of  existing  programs  written  by  the 
author,  with  modifications  being  included  where  necessary.  The  author 
makes  no  claim  to  originality  or  to  efficiency  of  operation  with  re¬ 
gard  to  the  program.  The  main  purpose  of  writing  such  a  program  was 
to  have  available  for  cost  analysts  an  easily  workable,  user-oriented, 
curve-fitting  computer  program--especially  adapted  to  handle  th  lath-  I 

ematical  functions  most  coranonly  used  in  the  development  of  estimating  | 

relationships.  j 
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SUMMA.RY 

This  Memorandum  describes  a  FORTRAN- IV  curve-fitting  compel  - 
progra  that  has  been  developed  within  the  kAND  Cost  Analysis  Depart¬ 
ment.  'flie  program  makes  least -squares  determinations  of  the  param¬ 
eters  of  any  of  five  mathematical  functions  selected  by  the  user, 
given  a  set  of  observations  on  the  dependent  and  independent  variables 
of  interest.  The  functions  available  in  the  program  are  the  line, 
parabola,  power,  asymptotic-power,  and  exponential.  Up  to  three  in¬ 
dependent  variables  may  be  used  for  the  line  and  power  functions. 

Also,  the  Y-intercept  may  be  specified  for  the  line,  parabola,  or 
asymptot ic -power  function. 

A  discussion  of  the  characteristics  of  the  functions  is  presented 
in  Section  I,  including  an  examination  of  those  nonlinear  functions 
that  require  special  methods  for  solution.  Also  included  is  a  brief 
discussion  of  the  statistics  used  in  the  program.  Specific  details 
on  the  operation  of  the  program  arc  presented  in  Section  II.  This 
section  also  treats  the  options  available  to  the  user.  Program  out¬ 
puts  are  discussed  in  Section  III.  For  the  benefit  of  the  reader, 
sample  outputs  from  two  ..ms  are  shown. 

Mathematical  considerations  relating  to  nonl inear-least -squares 
solutions  are  treated  in  Appendices  A  and  B.  A  listing  of  the  FORTRAN- 
IV  computer  progiam  is  presented  in  Appendix  C. 
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I.  INTRODUCTION 


PROGRAM  DESCRIPTION 

A  FORTRAN-IV  curve-fitting  computer  program  (CURVES)  has  been 
written  by  the  author  that  makes  least -squares  determinations  of  the 
parameters  of  any  of  five  types  of  functions,  given  a  set  of  observa¬ 
tions  on  the  dependent  and  independent  variables  of  interest.  These 
functions  are  commonly  used  in  the  derivation  of  cost  analysis  esti¬ 
mating  relationships,  and  arc:  (a)  line,  (b)  parabola,  (c)  power, 

(d)  asymptotic-power,  and  (e)  exponential.  They  are  described  in  de¬ 
tail  in  subsequent  parts  of  this  section.  Standard  statistics  re¬ 
lating  to  "goodness -of -f it"  measures  are  also  calculated  in  the  pro¬ 
gram.  No  predictive  statistics  are  included,  however,  because  of  the 
difficulty  of  obtaining  such  statistics  for  the  nonlinear  functions-- 
the  latter  three  above.  Consequently,  the  program  is  intended  essen¬ 
tially  for  curve-fitting. 

The  CURVES  program  can  handle  up  to  200  data  points  for  each  re¬ 
gression  and  is  so  structured  that  if  a  set  of  data  cards  contains 
data  for  several  separate  regressions,  that  set  needs  to  be  entered 
only  once.  This  obviates  the  need  for  duplicating  such  input  data 
decks  for  each  regression  run.  A  variable-format  procedure  is  provided 
the  user  so  that  data  may  be  entered  in  any  order  on  the  input  cards. 
Also,  an  option  is  provided  to  allow  the  user  to  specify  the  Y-intercept 
value  (regression  constant)  for  the  line,  parabola,  and  asymptotic- 
power  functions. 

The  program  is  written  completely  in  FORTRAN-IV,  using  A4  formats 
for  all  alphanumeric  information.  No  matrix-inversion  or  other  sub¬ 
programs  are  used.  All  solutions  are  made  through  either  standard, 
algebraic  methods  for  the  linear  and  parabolic  cases  or  through  iter¬ 
ative  methods  r:or  the  other  cases.  Cons equent ' y ,  the  program  should 
be  readily  adaptable  to  other  computer  systems. 

*$C 

In  this  Memorandum,  a  linear  function  is  defined  as  one  which 
is  linear  with  respect  to  all  of  its  parameters.  Under  this  definition, 
the  parabola  is  considered  to  be  a  linear  function. 
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FUNCT ION  TYPES 

The  functions  available  in  this  program  were  chosen  principally 
cn  the  basis  of  their  application  to  the  derivation  of  cost  analysis 
estimating  relationships.  They  are: 


ine  (containing  up  to  three  ind ependent  variables) 


Y  ’  A  +  B  •  (XI)  ,  \ 


Y  =  A  '  «  •  (XI) 

+  C  • 

CO, 

Y  =  A  +  B  •  (XI) 

+  c  • 

(X2)  +  L>  •  (X3> 

Y-intercep 

N  (A)  may  be 

spec  if ied . 

Parabola 

Y  -  A  +  B  ■  (XI) 

+  c  • 

2 

(XI)  .  j 

Power  (containing  up 

to  three  independent  variables) 

Y  *  A  *  (X1)B, 

Y  “  A  •  (X1)B  • 

(x?)c, 

Y  =  A  *  (X1)B  • 

<X2)C  • 

(X3)D. 

4 .  As ynp to tic- Power 


Y  =  \A  ■  (X1)B1  +  C 


Y-intercept  (C) 
may  be  specified. 


5 •  Sxponant la  1 

y  .  t,U  +  B  •  (XI)] 


•ff 

As  in  FORTRAN  notation  in  the  program,  the  independent  variables 
are  represented  by  XI,  X2 ,  and  X3,  respectively  When  only  one  inde¬ 
pendent  variable  Is  considered,  XI  is  used.  Also,  as  explained  later, 
only  positive  values  are  considered  for  all  variables . 


■ 


Y  -  dependent  variable, 

XI,  X2,  X3  *  independent  variables, 

A,  B,  C,  D  *  parameters  to  be  determined  by  least -squares  methods, 
e  =  constant  =.  2.71828. 


FUNCTION  CHARACTERISTICS 

Examples  of  some  of  the  types  of  curves  that  can  be  obtained  from 
the  five  functions  are  shewn  in  Fig.  i  for  a  one- independent -variable 


The  linear  form  is  the  most  simple  of  the  forms  treated  here.  Its 
characteristics  are  well  known  and,  in  the  opinion  of  the  author,  need 
no  further  elaboration.  Hie  user  has  the  option  of  using  up  to  three 
independent  variables  and  also  the  option  of  specifying  the  Y-intercept 
(A). 


Parabola 

Sometimes  the  parabolic  function  is  used  to  represent  points  that 
lie  along  a  curve  having  a  Y-intercept  (including  zero).  However,  one 
must  be  aware  that  since  this  function  is  actually  a  polynomial  of  de¬ 
gree  2,  its  curve  always  has  a  maximum  or  minimum  point  (vertex).  This 
means  that  the  effect  of  the  independent  variable  (XI)  on  the  dependent 
variable  (Y)  is  reversed  once  tuis  point  is  traversed.  Again,  the  user 
has  the  option  of  specifying  the  Y-intercept  (A). 


Power 

Hie  power  function  is  one  of  the  more  common  functions  used  in 
cost  analysis  work.  A  plot  of  its  logarithmic  counterpart,  the  log- 
linear  form,  is  known  as  the  "learning  curve"  or  "improvement 
cost  curve."  However,  for  reasons  discussed  later,  the  power,  rather 
than  the  logarithmic ,  form  is  used  in  this  program.  For  this  function, 


Y  A  ■  8  •  X I 


Parobo la 

Y  A  •  B  *  X I  *  C  •  X  r 


Power 
Y  A  •  X I ! 


Asymptotic -power 
Y  -  f  A  *  X !  8 )  *  C 


0  <  B  <  1 


Hort  zontal 
asymptote  (  Y  C 


Goes  through  origin  for  B>0 


Exponential 

Y  e  A  -f  B  *  X ' 


Fig  ! --Examples  oi  Curves  used  in  Program  for  a 
One-Independent -Variable  Case 
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the  user  has  the  option  of  using  up  to  three  independent  variables. 

For  positive  exponent  B,  this  curve  always  passes  through  the  origin, 
as  shown  in  the  figure.  Therefore,  it  should  never  be  used  where  a 
positive  Y-intercept  is  desired  cr  logically  required.  For  negative 
B,  the  curve  is  undefined  at  XI  =  0  and  is  a  declining  curve,  approach¬ 
ing  asymptotically  the  XI -axis  as  XI  becomes  large. 

Asymptot ic -Power 

An  examination  of  the  fourth  function,  the  asymptot ic -power ,  shows 
that  the  curve  has  a  horizontal  asymptote  of  Y  -  C  for  negative  B. 

g 

That  is,  as  XI  becomes  large,  the  first  term  (A  •  XI  )  approaches  zero, 
and  hence  the  value  of  Y  approaches  that  of  the  constant  term  C.  Con¬ 
sequently,  there  is  a  level-off  effect  associated  with  this  curve  for 
negative  B.  Tins  function  may  thus  be  used  to  represent  points  that 
lie  along  a  curve  either  increasing  or  decreasing  to  a  horizontal 
asymptote.  Like  the  power  curve,  this  curve  is  undefined  at  XI  =  0 
for  neg_.ive  B.  For  positive  B,  there  is  a  Y-intercept  equal  to  C. 

g 

As  XI  becomes  large,  the  first  term  (A  •  XI  )  ultimately  becomes  large 
compared  with  C,  and  therefore  Y  approaches  a  pure  power  function 
(A  •  X1B)  in  this  region  of  XI. 

A  plot  of  the  asymptot ic -power  function  on  log-log  paper  produces 
a  curved  line  at.  low  values  of  XI  that  approaches  cither  a  horizontal 
(B  <  0)  or  inclined  (B  >  0)  asymptote  at  high  values  of  XI.  Where  a 
positive  Y-intercept  is  desired,  the  user  may  specify  a  value  for  the 
constant  term  C. 

Exponent  fal 

The  last  form,  the  exponential,  is  used  to  represent  points  that 

A 

lie  along  a  curve  having  a  positive  Y-intercept  (e  ).  The  curve  may 
be  either  a  rising  (B  >  0)  or  falling  curve  (B  <  0) --the  falling  curve 
approaching  asymptotically  the  Xl-axis.  The  logarithmic  counterpart 
of  the  exponential  function  is  the  semi  log  function,  which  produces  a 
straight  line  on  semi  log  paper.  That  is.  In  Y  is  a  linear  function 
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‘V 

of  XI.  As  was  the  case  with  the  power  function,  for  reasons  discussed 
later,  the  exponential,  rather  than  the  semilogarithmic ,  form  is  re¬ 
tained  for  this  program. 


NONLINEAR-LEAST-SQUARES  SOLUTIONS 

It  can  be  shown  mathenetically  that  the  leas t -squares  solutions 
of  the  parameters  of  any  function  are  always  exact  and  unique  provided 
that  the  function  is  lineat  with  respect  to  all  of  its  parameters. 
Therefore,  for  this  program,  the  line  and  parabola  produce  exact  and 
unique  solutions.  (The  term  exact  is  used  to  refer  to  a  solution  that 
can  be  obtained  algebraically.)  However,  the  latter  three  functions 
are  not  all  linear  in  terns  of  their  parameters.  Thus,  their  solutions 
are  not  exact  and,  as  shown  later,  may  not  represent  absolute  minimums 
of  the  sum  of  squares  of  the  Y  residuals.  They  must  be  obtained  in 
some  other  way--usually  through  some  type  of  iterative  procedure.  (The 
general  principles  of  such  procedures  and  other  mathematical  consider¬ 
ations  relating  to  the  solutions  of  nonlinear-least -squares  equations 
are  presented  in  Appendix  A.) 

Foi  the  power  and  exponential  functions,  a  modified  Gauss-Newton 
method  is  used,  in  which  initial  estimates  are  obtained  from  the  log¬ 
arithmic  solutions  (which  are  exact)  and  then  correct iono ,  guaranteed 
to  produce  convergence  to  a  solution,  are  applied  to  those  initial  es¬ 
timates.  nils  procedure  is  repeated  until  the  absolute  change  in  the 

value  of  each  parameter  becomes  equal  to,  or  less  than,  some  predeter- 

.  -8  .  *■* 

mined  value  (10  in  the  program). 

The  solution  of  the  asymptotic-power  function  is  based  on  another 
type  of  iterative  procedu  because  there  appears  to  be  no  easy  way  to 


All  logarithms  discussed  herein  arc  natural  logarithms  (base  e) 
and  are  represented  by  !n. 

This  procedure  is  described  in  detail  in  KM-4879-PR  by  C.  A. 
Graver  and  H,  E.  Boren,  Jr.,  Multivariate  Logarithmic  £.:d  Exponential 
Regression  Models,  The  RAND  Corporation,  July  1967.  It  may  be  noted 
that  conceptually,  the  solutions  for  the  power  and  exponential  func¬ 
tions  are  each  different  than  for  their  logarithmic  counterparts  (see 
Appendix  A).  Also,  the  term  "exponential"  in  the  above  referenced  RM 
is  equivalent  to  it,  term  "power"  in  this  Memorandum. 


5 


-7- 


obtain  the  initial  guesses  that  are  required  for  the  modified  Gauss  - 
Newton  method,  nils  procedure  is  treated  in  Appendix  B. 

STATISTICAL  CONSIDERATIONS 

Because  of  the  difficulty  of  calculating  and  applying  predictive- 

type  statistics  for  the  nonlinear  functions,  it  was  decided  to  use  only 

ic 

"goodness  of  fit"  statistics  in  the  program.  Consequently,  this  pro¬ 
gram  should  be  regarded  as  essentially  a  curve-fitting  program  with 
only  those  statistics  being  used  that  relate  to  how  well  the  curve  fits 
that  particular  set  of  data.  Also,  it  should  be  noted  that  the  statis¬ 
tics  may  not  have  exactly  the  same  meaning  for  the  power,  asymptotic- 
power,  and  exponential  functions  as  for  the  line  and  parabola  because 
of  the  nonlinear  characteristics  of  the  former  three.  In  general, 
statistics  for  nonlinear  functions  should  be  used  with  care.  For  ex¬ 
ample,  unless  there  is  proof  to  the  contrary,  the  F  statistic  for  a 
nonlinear  function  probably  should  not  be  compared  with  the  F  table. 
Such  statistics  should  generally  be  used  only  qualitatively--not 
quant itat ively--unt il  a  thorough  investigation  is  made  into  the  appli¬ 
cation  of  such  statistics  to  nonlinear  functions. 

The  principal  reason  for  omitting  the  logarithmic  forms  in  this 
program  is  that  it  is  very  difficult  to  compare  their  fits  statisti¬ 
cally  with  those  of  the  nonlogar ithmic  forms  (see  Appendix  A).  As  a 
result,  no  logarithmic  curves  ore  used  and  the  statistical  resu’ts 
relating  to  the  five  functions  used  in  this  program  can  be  compared 
more  directly.  However,  since  the  iterative  solutions  for  the  power 
and  exponential  functions  require  that  their  logarithmic  solutions 
be  determined  for  the  initial  estimates,  these  solutions  are  also 
printed  in  the  output  (without  any  related  statistics)  for  the  benefit 
of  the  user. 

A  summary  of  the  statistical  equations  presented  in  Section 
III,  following  the  discussion  on  program  outputs. 

The  use  of  predictive  statistics  for  the  power  function  is  treat¬ 
ed  in  RM-4879  (see  previous  footnote).  However,  that  program  requires 
many  additional  subroutines,  which  in  the  opinion  of  the  author,  would 
make  the  CURVES  program  prohibitively  large,  slow  in  operation,  and 
restricted  to  50  or  less  data  points. 
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II.  INPUT  PROCEDURES 


The  flow  of  operations  within  the  program  is  depicted  in  Fig.  2. 
The  program  is  so  structured  that  many  sets  of  data  may  be  entered,  in 
which  each  set  200  data  points^  constitutes  a  run.  As  soon  each 
set  is  read  in.  the  program  operates  on  that  set  before  proceeding  to 
the  next  set  of  input  data.  Each  data  set  may  be  entered  on  a  separate 
deck  of  cards.  On  the  other  hand,  several  or  all  of  the  data  sets, 
if  space  on  the  cards  permits,  may  be  entered  on  one  deck  of  cards, 
thus  effecting  considerable  savings  in  the  use  of  cards  and  in  the 
effort  of  duplicating  a  deck  of  cards  containing  data  for  several 
runs.  A  variable  format  procedure  is  used,  allowing  much  flexibility 


in  the  format  of  the  input 

data . 

Listed  below  are  the 

types  of  cards  that  must  be  entered  for  the 

first  run. 

1. 

Title  card 

2. 

Order  card 

)  Need  to  be  entered  only  once  if  input  data 
/  for  all  runs  are  to  be  entered  in  same 

3. 

Format  card 

7  format 

4. 

Scale  card 

Used  only  if  data  are  to  be  scaled  or  if 
Y-intercept  is  to  be  specified. 

5. 

Data  cards 

6. 

Blank  card 

7. 

End  card 

Optional 

TITLE  CARD 

Hie  titl 

e  card  must  b' 

entered  for  each  run.  In  addition  to  the 

title  (alphan 

umerlc),  this 

card  also  contains  other  information  about 

the  run.  If 

the  two  cards 

relating  to  the  variable-format  procedures 

(order  and  format  cards)  are  to  be  read,  a  "1"  is  entered  in  Col.  1  of 
the  title  card.  For  th  first  run,  a  "1"  must  be  entered.  For  subse¬ 
quent  runs  involving  different  input  formats,  it  still  must  be  entered 
However,  if  data  for  all  subsequent  runs  are  to  be  entered  on  separate 
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START 


SET  PAGE  NUMBER  TO  1 
IPAGE  =  I 


SET  RUN  COUNTER  TO  1 
NRUN  -  1 


T 


PRINT  MAIN  HEADING^ 
f  AND  PAGE  NUMBER  ON 
FIRST  PAGE.  I 


1 


SUBROUTINE  INPUT 

READS  INPUT  DATA  FROM  CARDS 
OR  UTILITY  DISK  AND,  IF  SO 
DESIGNATED,  WRITES  DATA  ONTO 
UTILITY  DISK  FOR  REREADING 
DURING  THE  NEXT  RUN. 


SUBROUTINE  CHECK 

CHECKS  INPUT  DATA  FOR  ERRORS 
AND,  IF  SO  DESIGNATED, 
SCALES  DATA. 


STEP  PAGE  NUM8ER  BY  1 
IPAGE  =  IPAGE  +  1 


CLEAR  COMMON. 
(SET  VALUES  OF  MOST 
VARIABLES  IN  COMMON 
TO  ZERO) 


SET  SUBROUTINE 
INDICATOR  TO  1  . 
1ND  =  I 


SUBROUTINE  SUMS 

ORDERS  THE  DATA  FROM  LOW  TO 
HIGH  VALUES  OF  Y  IF  INDICATED; 

OBTAINS  VARIOUS  SUMS  AND 
MEANS  OF  THE  INPUT  DATA,  AND 
OBTAINS  THE  STANDARD 
DEVIATIONS  OF  THE  INPUT  DATA. 


©* 


Asy  mp  tot  i  c  -paw  e  r 


CHECK  VALUE 
OF  EQUATION 
DESIGNATOR, 
IEQ 


SUBROUTINE  READ 

1 

i 

SUBROUTINE  PRtNT 

READS  TITLE  CARD  A  0, 

PRINTS  APPROPRIATE 

IF  SO  DESIGNATED,  ORDER 

— 

\  SUBHEADING  FOR 

CARD,  VARIABLE-FORMAT 

/  EQUATION  BEING 

CARD,  AND  SCALE  CARD. 

USED. 

(Parabola) 


.fcolalQ 


*Eqaatior»  ured  in  program. 

I.  LINE 

Y  ♦  A  +  S  (XI  ) 

Y  =  A  ♦  B  (XI )  +  C  (X2) 

Y  =  A  ♦ B  (XI ) ♦ C  (  X2  )  •*  D  ( X3  ) 


2.  PARABOLA 

Y  =  A  *  B  (X!  )  ♦  C  (XI  )3 

3.  POWER 

Y  =  A  (XI)'  (X2)C  (X3)° 


Intercept 
A  may  be 
tpecified 


4.  ASYMPTOTIC -POWER 
Y  =[A  (XI)*]>C 

5.  EXPONENTIAL 

y  ( **)] 


Intercept 
C  may  be 
ipecified 


Fig.  2--F1ow  ol  Operations 
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Fig.  2--Flo%  of  Operations  (Cent.) 
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Fig.  2 --Flow  of  Opera tlcms  (Cent.) 
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docks  of  cards  using  iho  same  format  as  the  first  run,  then  the  "1" 
only  needs  to  be  entered  for  the  first  run.  Whenever  the  1  is  not 
entered,  then  the  order  and  format  cards  are  not  entered. 

If  input  data  representing  several  runs  are  entered  on  one  deck 
of  cards  so  that  for  subsequent  runs  those  same  cards  will,  in  effect, 
be  reused,  and  hence  re-read  (in  different  fields),  then  a  "1"  is  enter - 
ed  in  Col ■  2  of  the  title  card  for  the  first  set  of  such  data  to  be 
read  in.  This  causes  the  machine  not  only  to  read  the  first  set  of 
data  but  to  write  all  of  the  input  data  onto  a  utility  disk  for  re¬ 
reading  during  subsequent  runs.  Unless  the  user  selects  another  util¬ 
ity  disk  or  tape,  the  program  automatically  uses  utility  disk  S.SU04 
(FORTRAN  logical  unit  4)  for  this  operation. 

For  the  remaining  runs  that  use  the  input  data  from  the  same  deck 
of  cards,  a  "2"  is  entered  in  Col ■  2  of  each  title  card  for  those  runs. 
This  causes  the  machine  to  read  the  input  data  from  the  utility  disk 
(instead  of  from  new  cards)  in  accordance  with  the  format  instructions 
so  entered. 

Column  3  is  used  for  the  function  designator.  An  integer  from 
"1"  through  "5"  is  entered  to  designate  which  function  is  being  con¬ 
sidered  for  that  run.  The  integer  designators  are  as  follows: 

Integer  Function 

1  Line 

2  Parabola 

3  Power 

4  Asymptotic-power 

5  Exponential 

There  must  be  an  integer  of  one  of  the  above  values  entered  in  Col. 

3  for  the  first  run.  If  Col.  3  is  left  blank  after  the  first  run, 
then  the  value  for  the  previous  run  is  used.  Thus,  if  the  same  type 
of  function  is  being  examined  for  a  series  of  runs,  its  designator 
needs  to  be  entered  only  for  the  first  run. 

Column  4  is  used  to  designate  whether  the  input  data  are  to  be 
ordered  from  low  to  high  values  of  Y.  A  value  of  "1''  signifies  that 
the  data  are  to  be  ordered.  For  the  first  run  a  blank  (zero)  signifies 
that  the  data  are  not  to  be  ordered.  However,  for  subsequent  runs  a 
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blunk  (zero)  signifies  to  the  machine  that  the  value  of  the  order 
designator  for  the  preceding  run  is  to  be  used.  Again,  this  is  done 
so  that  if  all  runs  in  a  scries  are  to  be  either  ordered  or  unordered, 
the  order  designator  will  only  have  to  be  entered  for  the  first  run. 

In  the  case  where  a  zero  is  desired  fur  the  designator  after  a  "1" 
has  been  entered  previously,  a  "2"  must  be  entered.  This  in  effect 
sets  the  value  of  the  designator  to  zero. 

The  scale  designator  is  entered  in  Col.  5.  If  any  of  the  data 
are  to  be  scaled  (using  the  Scale  Card  described  later),  a  "1"  is 
entered  in  this  column.  Otherwise  it  is  left  blank  (zero). 

Column  6  is  used  to  designate  whether  the  Y-intercept  term  for 
the  line,  parabola,  or  asymptotic-power  function  is  to  be  specified. 
This  is  done  by  entering  a  "1."  Otherwise  Col.  1  is  left  blank  (zero). 
If  a  "l”  is  entered  in  either  Cols.  5  or  6,  or  both,  then  the  scale 
card  is  entered  in  the  order  shown  previously. 

Columns  9  through  72  are  reserved  for  the  title.  The  title  may 
consist  of  any  alphanumeric  symbols. 

A  summary  of  the  information  on  the  title  card  is  given  in  Table 
1.  An  example  of  a  title  card  is  shown  in  Fig.  3,  in  which  a  linear 
regression  is  to  be  made  on  the  input  data  ("1"  in  Col.  3).  The  data 
,ire  to  be  ordered  with  respect  to  Y  ("1"  in  Col.  4)  and  are  to  be 
scaled  ("1"  in  Col.  5).  The  "1"  in  Col.  1  indicates  that  the  order 
and  format  cards  are  to  be  read  next. 

ORDER  CARD 

The  next  card  (when  used  as  specified  by  a  "1"  in  Col.  1  of  the 
title  card)  indicates  the  order  (from  left  to  right)  in  which  the  data 
are  located  on  the  data  cards.  The  order  (designated  by  alphanumeric 
symbols)  is  entered  in  Cols.  1-2,  4-5,  7-8,  10-11,  and  13-14.  Depend¬ 
ing  on  the  number  of  independent  variables  being  used  and  on  whether 
identifiers  are  being  used,  Cols.  7-14  may  not  be  required.  The  sym¬ 
bols  used  to  show  the  order  are  as  follows: 
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i 

Table  1 

SIMMARY  OF  INFORMATION  ON  TITLE  CARD 


Columns 

Remarks 

1 

A  "1"  indicates  that  the  order  card  and  format 
card  are  to  be  read,  respectively,  following  the 
title  card  If  blank,  no  such  cards  are  to  be 

read.  A  "1"  roust  be  entered  for  the  first  run. 

2 

A  "1"  indicates  that  the  input  data  cards  contain 
data  for  several  runs  and  are  to  be  written  onto  a 
utility  J'sW  (Sysr  n  Unit  S.Sl'04,  FORTRAN  Logical. 
Unit  4),  A  *'2 ”  indicates  that  the  input  data  are 
to  be  read  from  the  disk.  If  blank,  the  input 
data  cards  for  this  run  are  to  be  read  only  once- 

3 

II 

An  integer  from  "1"  to  "3"  is  used  to  designate 
which  function  is  being  used.  This  is  done  as 
follows . 

1  --  Linear 

2  --  Farabolic 

3  --  Fower 

4  --  Asymptot ic -power 

5  --  Exponential 

If  blank,  after  the  first  run,  the  value  for  the 
preceding  run  is  used. 

4 

A  "1"  Indicates  that  the  data  are  to  be  ordered 
from  low  to  high  values  of  Y.  For  the  first  run 
a  blank  Indicates  that  the  data  are  not  to  be 
ordered.  If  blank  for  subsequent  runs,  the  value 
for  the  preceding  run  is  used.  A  "2"  is  used  for 
subsequent  runs  to  restore  the  designator  to  zero 
when  desired. 

5 

A  "1"  indicates  that  the  data  are  to  be  scaled. 
Otherwise,  it  is  left  blank. 

A  "1"  indicates  that  a  Y-inte^  .  pt  is  to  be 
specified  for  either  the  linear,  parabolic,  or 
asymptot ic -power  case.  Otherwise,  it  is  left 
b lank . 


6 
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Table  1  (Coni.) 

SIMMARY  OF  INFORMATION  ON  TITLE  CARD 


9-72 


Title  for  run.  May  consist  of  any  alphanumeric 
symbols . 
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Fig.  3--ExatnpIr>  of  TitK  Card 
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Symbo 1 

ID  . , 

Y1  . 

XI  . 

X2  . 

X3  . 


Type  of  Data 

Identifier  (a Iphanutnerir)  (optional) 
Dependent  variable  (required) 

First  independent  variable  (required) 
Second  independent  variable  (optional) 
Third  independent  vriable  (optional) 


Hay  be  in  any 
order  from  left 
to  right. 


Suppose  that  a  sp*  of  data 
the  three  independent  variables 
in  Cols.  1-12,  13-24,  25-36,  and 
that  an  identifier  (a  six-digit 
"XI,"  "X2 ,  "  •’X3,"  "Y  2,"  and  "ID" 
Cols.  1-2,  4-5,  7-8,  10-11,  and 


is  to  be  entered  in  which  values  for 
and  the  dependent  variable  are  located 
37-48,  respectively.  Suppose  also 
integer)  is  in  Cols.  55-60.  Then, 
would  be  entered  respectively,  in 
13-14  to  show  the  above  order  across 


the  card. 


All  alphanumeric  information  is  treated  in  A4  formats  in  this 
program  in  order  to  be  adaptable  to  IBM-360  systems.  In  addition,  the 
identifiers  may  be  entered  in  either  A4  or  2A4  formats.  This  is  in¬ 
dicated  in  Col.  16  on  the  order  card  by  either  a  "1"  (A4)  or  a  "2" 
(2a4),  If  Col.  16  is  left  blank,  after  the  first  run,  then  the  value 
for  the  preceding  run  is  used. 

The  order  card  is  shown  in  Fig.  4  for  the  above  example,  in  which 
a  2A4  format  is  to  be  read  for  the  six-digit  identifier.  Since  the 
commas  separating  the  symbols  are  in  columns  that  are  not  read  by  the 
machine,  they  may  be  used  for  the  purpose  of  clarification. 


FORMAT  CARD 

rhe  format  card  indicates  where  the  data  are  located  on  the  data 
cards.  Again,  this  card  is  used  only  if  a  "1"  is  entered  in  Col .  1 
of  the  title  card.  This  card  must  begin  with  a  left-hand  parenthesis 
and  end  with  a  right-hand  parenthesis  before  ,  and  the  infor¬ 

mation  within  the  parentheses  must  conform  to  the  rules  for  FORTRAN 
formats.  Except  for  the  alphanumeric  identifiers,  all  input  data 
roust  be  in  real-number  (floating  point)  formats.  The  format  card 
shown  in  Fig.  5  would  be  used  for  the  previous  example. 
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Fig.  4- -Example  of  Order  Card 
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Fig.  5--Example  of  Format  Card 
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!  SCALE  CARD 

L  — - 

As  designated  by  a  "1"  in  either  Cols.  5  or  6  (or  both)  of  the 
title  card,  the  scale  card  is  used  either  if  the  data  are  to  be  scaled, 
or  if  the  Y~inter<_ept  is  to  be  specified  for  the  linear,  parabolic, 
or  asymptotic-power  function. 

Data  Scaling 

The  first  four  sets  of  two  columns  each  on  the  scale  card  are 
used  for  scaling  the  input  dependent  and  independent  variable  data 
when  required*  These  scale  indicators  (integers)  are  located  as  .'ol¬ 
io  vs* 


Column  Location  of 

Variable  Scaled  Scale  Designator 

Y  (dependent  variable)  1-2 

XI  (first  independent  variable)  3-4 

X2  (second  independent  variable)  5-6 

X3  (third  independent  variable)  7-8 

An  integer  (fixed-point)  number  is  used  to  indicate  the  number  of 
places  that  the  decimal  is  to  be  moved.  A  positive  number  indicates 
that  the  decimal  is  to  be  moved  to  the  right  as  many  places  as  is  the 
value  of  the  number.  A  negative  number  indicates  that  the  decimal  is 
to  be  moved  to  the  left  as  many  places  as  is  the  absolute  value  of  the 
number.  For  example,  suppose  r^ar  a  "3"  is  entered  in  Col.  2  (must  be 
right-justified).  Then  each  input  Y  value  will  have  its  decimal  point 
moved  three  places  to  he  right,  e.g.,  50.1  23  — ►50123.,  0.6127— -► 

612.7,  etc. 

Care  must  be  taken  to  enter  all  positive  integers  in  the  right- 
hand  column  of  the  two-column  set.  If,  for  example,  a  "3"  were  enter¬ 
ed  in  Col.  1  instead  of  Col.  2,  the  machine  would  read  the  number  as 
"30"  Instead  of  "3."  Any  scale  factor  entered  applies,  of  course,  to 

the  entire  set  of  the  corresponding  variable  for  that  run  only.  An  ' 

example  of  this  card  is  shown  in  Fig.  6,  in  which  the  first-independent-  I 

variable  data  (XI)  arc  «_o  be  scaled  down  by  a  factor  of  10.  j 

i 

* 

Data  may  also  be  scaled  by  using  P-fonnats. 

] 
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Fig.  6--Example  of  Scale  Card 
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Speclfled  Y-Intercept 

The  specified  Y-intercept  (regression  constant)  is  entered  in 
Cols.  11-20,  It  is  entered  in  real-number  form  anywhere  within  these 
columns.  The  implied  decimal  point  location  is  at  the  right  end  of 
the  field,  between  Cols.  20  and  21. 

DATA  CARDS 

Each  data  card  must  contain  at  least  a  pair  of  values--one  for 
the  dependent  variable  (Y)  and  one  for  the  independent  variable  (Xl). 
Each  data  set  constituting  a  run  must  contain  one  more  data  point  than 
the  number  of  parameters  being  solved  and  may  contain  up  tc  200  data 
points.  The  location  of  the  data  on  the  cards  must  be  in  exact  agree¬ 
ment  with  the  information  entered  on  t lie  order  and  format  cards,  or 
else  the  data  will  not  be  read  properly.  The  numerical  data  (dependent 
and  independent  variable  values)  must  be  entered  as  real  numbers,  and 
the  identifiers  (if  used)  as  alphanumeric  data. 

If,  for  any  data  card,  either  the  Y  field  or  any  of  the  X  fields, 
that  are  read,  but  r.ot  all,  are  blank  or  contain  zeros,  that  card  is 
skipped.  However,  if  all  X-Y  fields  are  blank  (zeros),  the  reading 
of  input  data  for  the  run  is  terminated  at  that  point  (see  Blank  Card 
below) . 

In  conjunction  with  the  above,  zero  or  negative  values  ate  not 
allowed  for  the  X-Y  input  data.  Another  reason  for  this  is  that  for 
the  power  and  as ymptot ic -power  functions  involving  negative  fractional 
exponents  of  negative  X  values,  such  expressions  are  meaningless. 

Also,  logarithms  are  used  in  the  solutions  of  the  nonlinear  functions, 
and  the  use  of  logarithms  restrict  the  input  data  to  values  greater 
than  zero. 

An  example  of  a  data  card  containing  data  in  the  format  depicted 
in  Figs,  4  and  5  is  shown  in  Fig  7. 

BLANK  CARD 


Each  pack  of  data  cards  must  always  end  with  a  blank  card.  This 
card  is  used  to  terminate  the  reading  of  the  input  data  for  a  given  run. 
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Fig.  7 -’Example  of  Data  Card 


END  card  (OPTIONAL) 


The  program  Is  so  structured  that  after  each  set  of  data  is  read 
In  and  processed,  the  machine  attempts  to  read  in  another  set  of  data 
After  the  final  run  when  the  machine  attempts  to  read  in  another  set 
of  data  that  does  not  exist,  control  is  returned  to  the  FORTRAN  monitor, 
ending  the  series  of  runs.  At  this  point,  a  statement  is  printed  as 
fol lows : 


End-of-data  encountered  on  svstem  input  file. 

In  order  to  provide  a  positive  method  of  terminating  a  series  of 
runs,  this  optional  method  is  provided  trie  user.  If  a  "9"  is  entered 
in  Col.  1  on  a  card  immediately  following  the  blank  card  of  the  last 
set  of  data  (actually  the  next  title  card),  it  will  c«  se  the  program 
to  Call  Exit  and  terminate  the  series  of  runs  at  that  ”>o;nt.  If  this 
is  done,  the  End-of-Data  statement  will  not  he  printed. 

SUMMARY 

Figure  8  shows  the  order  of  two  data  decks  for  a  series  of  three 
runs.  In  this  figure,  the  second  deck  of  data  also  contains  data  for 
the  third  mu.  Hence,  it  is  to  he  written  onto  a  utility  disk  for  re¬ 
reading  during  the  third  run. 


Third  run  i  Data  for  this  run 
are  to  be  read  off  disk  con¬ 
taining  all  data  from  second 
deck  of  cards.) 


Second  run  (  Data 
for  this  run  are  to 
be  read  in  same 
order  and  format  as 
first  run .  In  addition, 
the  cards  are  to  be 
written  c--.ro  ais!> 
since  they  also 
contain  data  for 
the  third  run  .  1 


Format  card 


Order  card 


Title  card  for 
third  run 


Blank  card  at 
end  of  second 
data  deck 


/  Data  deck  for 
second  and 
third  rut's 

f  Title  card  for 
1  second  run 

8lank  card  at  end  ~ 

of  first  data  deck  ) 


Data 

a  run 


points  for 
<  200) 


Data  deck  for 
first  run 


Format  carl 


Ordc  card 


/ 

Title  ca 

d  for 

first 

v'l 

tr>f  run 


k atn-  ru  of  Two  Da t a  Card  Packs  tor  Hire.-  Runs 


-25- 


III.  PROGRAM  OUTPUTS 


Examples  of  program  outputs  for  the  line  and  power  functions 
are  shown  in  Figs.  9-12.  The  particular  data  shown  have  no  meaning 
whatsoever  and  are  presented  only  for  the  purpose  of  displaying  the 
outputs.  Since  the  headings  are  rather  self-explanatory,  and  the 
statistics  have  already  been  discussed  in  Chapter  I,  no  explanation 
of  them  will  be  given.  It  should  be  emphasized  again  that  the  statis¬ 
tics  presented  may  have  different  meanings  with  regard  to  the  nonlinear 
functions  than  with  the  linear  functions.  They  should  not  be  used  as 
predictors  and  should  be  regarded  only  as  indicators  of  the  "goodness 
of  fit"  for  that  particular  data  set. 

As  is  shown  in  Appendix  A  for  the  power  and  exponential  functions, 
if  there  is  a  solution  of  the  parameters  lying  in  the  region  defined 
by 


2 


-  Y.  >  0, 

l 


where  Y  =  computed  value  of  Y  at  each  point  i, 

Y^  =  observed  value  of  Y  at  each  point  i, 

then  that  solution  represents  an  absolute  minimum  of  the  sum  of  squares 
of  Y  residuals  and  is  the  only  solution  in  that  region.  In  addition, 
it  seems  reasonable  to  assume  that  the  above  condition  should  hold  for 
those  functions  in  order  to  have  a  "good"  fit.  If  the  above  condition 
for  the  region  does  not  hold,  a  message  is  printed  noting  that  fact. 

For  the  parabolic  function,  the  X-Y  coordinates  of  the  maximum 
or  minimum  point  (vertex)  are  printed.  After  this  point  is  passed, 
the  effect  of  XI  on  Y  is  reversed. 

Figure  13  is  a  listing  of  the  input  data  for  the  two  runs  whose 
outputs  are  shown  in  Figs.  9-12.  The  data  are  in  the  same  format  on 
the  cards  as  shown  in  Figs.  4-7.  As  Fig.  13  shows,  the  input  data 
for  the  second  run  (power  function)  are  read  in  the  same  format  as 
for  the  first  run.  Therefore,  the  order  and  format  cards  are  not 
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85.28716 

27.  182  79 

71.16251 

100.  1  02  73 

145.27168 

1 12. 1H882 

14  7.2355  7 

25. 18892 

38.29910 

2  7.  1  802  7 

75  .  2  7 ' 64 

09.26615 

5.27715 

45.23519 

1  00  "•  04  ^ 

Title  card 
Order  card 
Format  card 
Scale  card 

27.94610 

00.09592 

101313 

1U0. 47230 

1 0  2  1  2  1 

09.9267*. 

101717 

33.60992 

100505 

03.40204 

101414 

40. 71304 

100797 

94.84002 

101515 

54.53670 

1 0  1  0  1  0 

107. 75160 

102  32  3 

4C.9791 7 

1  0(i  81!  8 

19.44080 

1  CO  202 

l  37. C 

102  72  7 

88. 2662  7 

101616 

24.06789 

ICO  10  3 

^  Input  data 

1  4  7.550  76 

1  0  30  1C 

14. 39608 

1  C  0  1  0  1 

111.  3524  7 

102424 

92.51133 

101919 

51.71629 

10090  ) 

40.40560 

100606 

94. 38172 

102020 

104.27262 

102222 

7  1.0968  7 

101212 

1  31.24108 

102626 

9C. 28862 

101818 

59.0 

1 0  1  1 1  l 

l  38.4614  7 

102828 

122.6934  7 

102525 

1 4  3.  72060 

10292 

Blank  card 

15  1. 20949 

200606 

1510.21359 

702929 

66  1.  12500 

201515 

880.20145 

2  02  0  2  0 

535.55565 

201313 

8C0.0 

201018 

60.25169 

200101 

530. 19078 

201111 

900. 76950 

202222 

1138.31207 

202625 

972.57176 

202121 

115.19201 

200702 

460. 16664 

200909 

1  14  5.  540  35 

202727 

1159.45163 

202626 

14  3.  0400  7 

2  00  30  3 

1  70. 90614 

2  00  70  7 

14  72.22902 

2u?82« 

7  79.  62800 

201  71  7 

425.47155 

2  c  0  h  n  0 

342. 1 5482 

2011505 

loco. 58891 

2  02  32  3 

46  1.9  784  7 

20  l  0 ' 0 

0  1  4.0  790 e 

201 919 

l  1  25.6763  7 

202424 

755. 79487 

201616 

5  31.  15436 

2nl 212 

265C. 01076 

2  0  3  0  10 

2  0  7.960  3  7 

7 C 0  *.04 

58*  70353 

701414 

Blank  card 

Fig.  13--Li9ting  of  Inputs  for  Two  Runs  whose 
Outputs  are  Shown  in  Figs.  9-12 
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required  fur  that  second  run.  Also,  there  is  no  scale  card  for  the 
second  run  because  none  of  the  data  are  scaled. 

A  summary  of  all  of  the  statistics  that  are  used  in  the  CURVES 
program  is  presented  in  Table  2  in  simplified  form.  One  should  be 
aware  that  when  the  Y-fntercept  is  specified,  not  only  are  the  degrees 
of  freedom  changed,  but  the  curve  ls  no  longer  forced  through  the 
means  of  the  observed  data.  Consequently,  the  statistics  will  have 
different  values.  In  addition,  for  the  linear  case  with  the  speci¬ 
fied  intercept,  the  F  ratio,  which  is  based  on  the  distribution  of  the 

★ 

Y  residuals,  may  have  its  distribution  altered.  Therefore,  a  compar¬ 
ison  of  its  value  with  tabulated  values  based  on  an  assumed  standard 
F  distribution  may  not  always  be  valid.  As  a  result,  F  ratios  for  the 
specified  Y-intercept  cases  may  not  always  be  comparable  with  those 
for  unspecified  Y-intercept  cases.  A  statement  to  this  effect  ls 
printed  in  the  Summary  Table  of  the  output  whenever  the  Y-intercept 
is  specified. 

A  listing  of  the  FORTRAN-IV  computer  program  Is  presented  in  Ap¬ 
pendix  C. 


As  was  noted  previously  in  Section  I  under  Statistical  Consider¬ 
ations  ,  F  ratios  for  nonlinear  functions  may  not  be  valid  for  compari¬ 
sons  with  tabulated  values. 


Table  2 

STATISTICAL  EQUATIONS  USED  IN  PROGRAM8 


Statist ic 


Sum  of  Squares  of  Y 
Deviations 


Equation 


N 

S1  -  Z  <Yi  -  YC1>: 


Standard  Error  of  the  Es¬ 
timate  of  Y  (Adjusted) 


Coefficient  of  Variation 
(Percent) 


Coefficient  of  Determination 
(Unadjusted) 


where 

N  ■  number  of  data  points 

Y^  ■  observed  Y  value  at  each 
point  i 

Y  *  calculated  Y  value  at  each 
C  point  i 


1  N-P 


where 

P  3  number  of  parameters  to  be 
solve  A 


ft) 


where 

Y  a  mean  of  observed  Y  values 


S .  «  1 - 

U  N 


)  Of.  -Y)2 

L  ...  1 


Table  2  (Cont . ) 

STATISTICAL  EQUATIONS  USED  IN  PROGRAM 


Statistic 


Mean  of  Absolute  Percent 
Y  Deviations 


Standard  Deviation  of 
Input  Variables 


Equation 


<vt  -v)2 


where 


V  *  value  of  any  input  vari- 
1  able 

V  *  mean  value  of  input  vari¬ 
able  V 
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Appendix  A 

NONLINEAR -LEAST-SQUARES  CONSIDERATIONS 

LOGARITHMIC  AND  NONLOGARITHM  I C  SOLUTIONS 

The  usual  procedure  for  making  least -squares  determinations  of 
the  parameters  of  the  power  or  exponential  function  is  to  first  con¬ 
vert  the  function  into  a  log-line.tr  (or  semi  log -linear)  function  One 
then  has  a  logarithmic  linear  equation  which  can  be  solved  directly 
through  use  of  the  standard  normal  least -squares  equations  for  t bi¬ 
linear  case.  However,  one  should  be  aware  that  such  a  solution  is  not 
the  same  conceptually  as  the  least -squares  solution  of  the  function 
before  it  is  converted  into  its  logarithmic  counterpart.  This  may  in 
seen  by  considering,  for  examp.o,  the  power  function  and  its  icuur'thr  c 
form . 

Let 

Y  -  A  •  (X1)B  ■  (X2)C  •  (X3)°, 

and 

In  Y  -  In  A  +  B  ■  In  XI  +  C  •  in  X2  +  D  •  In  X3 

For  a  least-squares  solution,  one  is  interested  in  minimizing  the 

* 

sum  ot  squares  of  the  Y  residuals  (denoted  by  Q).  Therefore,  for  the 
power  function, 

N 

h'  “  (Y1  -  Ycl)2  -  min.  , 

i-1 

and  for  the  logarithmic  function 

N 

Q’  *  (In  Y.  -  In  Y  )2  *  min.  , 
i  ci 

i-1 


* 

Throughout  this  discus; ion,  Q  is  used  to  represent  the  sum  of 
squares  of  the  Y  residuals. 
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or 


Q' 


°  min.  , 


where  N  «  number  of  data  points, 

■  observed  value  of  dependent  variable  at  each  point  i, 
Yci  ■  calculated  value  of  dependent  variable  at  each  point  i. 


In  the  logarithmic  case,  the  sum  of  squares  of  the  actual  dif¬ 
ferences  (residuals)  between  the  observed  and  calculated  Y  values  is 
not  being  minimized — rather  the  sum  of  squares  of  the  logarithms  of 
the  ratios  of  those  values  is  being  minimized.  Depending  on  the  ob¬ 
served  data,  Q  and  Q'  may  produce  significantly  different  solutions 
for  the  parameters  A,  B,  C,  D. 

It  may  also  be  seen  that  any  statistic  based  on  the  sum  of  squares 
of  Y  residuals,  such  as  the  coefficient  of  correlation,  may  be  mis¬ 
leading  if  used  to  compare  the  logarithmic  form  with  its  nonlogarithmic 
counterpart.  For  the  logarithmic  form,  such  statistics  are  in  loga¬ 
rithms  end  hence  have  different  meanings. 

The  question  as  to  whether,  say,  the  power  function  or  its  log¬ 
arithmic  form  is  more  appropriate  for  a  set  of  data  is  beyond  the 
scope  of  this  Memorandum.  An  answer  to  such  a  question  depends  on 
many  factors  including  the  errors  associated  with  the  data  and  what 
criterion  is  used  for  a  "good  fit."  For  the  interested  reader,  this 
question  is  treated  in  more  detail  in  RM-4879-PR. 


NON-LINEAR  SOLUTIONS 

It  is  a  necessary  condition  that  the  first  partial  derivatives 
of  Q  with  respect  to  the  parameters  must  be  zero  in  order  that  Q  be 
minimized.  This  is  not,  unfortunately,  a  sufficient  condition  for  a 
function  that  is  not  linear  with  respect  to  all  of  its  parameters. 

The  reason  for  this  is  that  if  Q  could  be  graphed  (in  c  llti-dimensional 
space)  for  this  type  of  function,  there  might  be  other  critical 

Graver,  Boren,  o£.  cit . 
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points --such  as  saddle  points  or  relative  maxima  or  minima  points -- 
where  the  first  partial  derivatives  would  also  be  zero.  A  test  that 
checks  for  this  possibility  is  the  examination  of  the  matrix  of  second 
partial  derivatives  of  Q,  which  is  a  generalization  of  the  second- 
derivative  test  for  a  one-parameter  case.  If  this  matrix  is  positive- 
definite  for  all  parameters  in  a  region  containing  a  solution,  it  can 
be  shown  that  the  solution  represents  an  absolute  minimum  for  Q  in 
that  region  and  is  the  only  solution  in  that  region.  However,  if  the 
matrix  is  not  positive-definite  in  that  region,  then  there  may  be 
other  "solutions"  for  the  same  set  of  data. 

With  regard  to  the  power  and  exponential  functions,  it  can  be 
shown  that  the  matrix  of  second -partial  derivatives  of  Q  will  be 
positive-definite  for  values  of  the  parameters  which  lie  in  the  re- 
gion  defined  by 

2  •  Yci  -  >  0,  (coefficient  A  >  0  for  power  function) 

where,  again,  Yci  *  calculated  value  of  dependent  variable  at  point  i, 
Y^  “  observed  value  of  dependent  variable  at  point  i. 

If  there  is  a  solution  in  the  region  defined  above,  then  that  solution 
represents  an  absolute  minimum  of  Q  and  is  the  only  solution  in  that 
region.  In  addition,  it  seems  reasonable  to  assume  that  the  above 
condition  should  hold  in  order  to  have  a  good  fit  for  the  power  and 
exponential  functions. 

In  summary,  one  should  be  aware  that  for  a  nonlinear  function 
the  "solution"  obtained  may  not  represent  an  absolute  minimum  for  Q. 

The  only  sure  way  to  know  is  to  try  all  combinations  of  the  parameters 
for  each  data  sample  to  determine  all  "solutions"  and  to  then  determine 

*9c 

H.  0.  Hartley,  "The  Modified  Gauss -Newton  Method  for  the  Fitting 
of  Non-linear  Regression  Functions  by  Leas t -squares Technometrics. 
Vol.  3,  No.  2,  May  1961,  pp.  273-274. 

fcic 

The  proof  is  given  in  Appendix  D  of  RM-4879-PR  by  Graver,  Boren, 
op.  cit . 
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which  sol ut ion  gives  tin-  1  west  sum  of  squares  of  Y  residuals .  for 
practical  reasons  this  is  very  difficult  to  do.  However,  otn  must 
rt  mber  that  he  is  attempting  to  find  a  solution  to  a  function  that 
adequately  represents  the  data.  Whether  or  not  there  are  solutions 
in  other  unknown  regions  may  be  rather  unimportant  if  the  solution 
that  is  found  is  satisfactory  to  the  analyst --that  is,  it  satisfies 

'k 

his  criterion  for  a  go o'  fit. 

■k 

For  further  information  on  ^  'nl  inear  least -squares  solutions 
the  reader  is  referred  to  Applied  Regression  Analysis,  by  N.  R.  Draper 
rod  11.  Smith,  John  Wiley  and  Sons,  Inc.,  New  York,  London,  Sydney, 
196b,  Chap.  10,  pp.  263-304. 
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Appendix  B 

LEAST-SQUARES  SOLUTION  FOR  AS YMPTOTIC- PiWER  FUNCTION 


DETERMINATION  OF  PARAMETERS  A,  B,  C 


For  this  program,  a  least -squares  fit  is  assumed  for  the  data, 
in  which  the  sun  of  the  squares  of  i.he  Y  residuals  (differences  between 
the  observed  values  of  Y  ar.d  the  corresponding  calculated  values  of 
Y  from  the  regression  equation)  is  a  minimum.  Therefore,  the  Y  resid¬ 
ual  (R)  at  any  point  i  is 


or 


Ri  -  Yi  -  (A  •  XI*  +  C), 


(1) 


where 


ci 

X1i 
A,  B,  C 


Y  residual  at  point  i, 
observed  value  of  Y  at  point  i, 
calculated  value  of  Y  at  point  i, 
value  of  independent  variable  at  point  i, 
parameters  to  be  determined. 


The  requirement  for  a  le *st -squares  fit  for  N  sets  of  points  is  that 
the  sum  of  the  squares  of  the  Y  residuals  (denoted  by  Q)  shall  be  a 
minimum;  that  is, 

N 

Q  *  y  (Yt  -  A  •  XI®  -  C)2  -  min.  (2) 

i-1 

g 

First  let  Zj  ■  Xl^.  Then  Equation  (2)  can  be  written  as  follows: 
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i 

! 

I 


N 

Q  =-  (Yt  -  A  *  Zi  -  C)2  =  min., 
i=l 

* 

From  normal  regression  equations 

N 


(Y.  -  Y)  •  (Zt  -  Z) 


(Z .  -  Z)2 


i=l 

C  =  Y  -  A  •  1  . 


(3) 


(4) 


(5) 


Substituting  the  above  expression  for  C  into  Equation  (3)  gives: 


N 

Q  =  V  [Yt  -  A  •  Zt  -  (Y  -  A  •  Z)]2  , 
i-1 

N 

Q  -  2  r<Yt  -  Y)  -  A  •  <ZA  -  z)]2  , 

1=1 

N  N 

Q  =  /  <Y.  -  Y)2  -  2  •  A  •  7  [(Y.  -  Y> 

L-->  J*  1 

i*i  i=i 

N 

.  (Zt  -  Z)]+  A2  •  Y  (Zi  -  Z)2  # 

i=l 


S  -  V  (Y  -  Y)2  , 
yy  L,  i 

1=1 


Overscores  are  used  to  denote  means. 


N 


S 

yz 


L  (Yi  “  Y>  *  <zi  -  z> 

i-i 


N 


1-1 


Therefore, 


Q 


S 

yy 


s 

zz 


fr<w  Equation  (4)  and  the  above  definitions; 


S 


zz 


Thus , 


Q  -  S  -  2 

yy 


•  s 

zz 


9 


(6) 


However , 


Q 


s 

yy 


2 

(s  y 
y* 
s 

zz 


N 


S 

yz 


i-1 


i  =  l 


-  Y) 


(7) 


i 


I 

1 


-  Z  • 
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N 


N 


3  =  (Y.  -  Y)  -  Z.,  since  (Y.  -  Y)  =  0. 

yz  l  7  1  .  -  l 

i=l  i=l 


Also, 


s«  •  (zi  -  z>  • 


i-i 


-  2  .  Z 


zz 


i=l 


N  JN 

z.  +  z2, 

l 

i=l  i=l 


zz 


Zt  -  2 
1 


1 

1 


(  N 


4  N  '  (n”  *  I  ?'i 

V  i-1  ) 


N 

.1  Zi 

i=l 

1 2 


•  /  z- 

L  i 

i-1 


N 

-  V  72  1 

'zz  ‘  i  '  N 

i=l 


Finally,  using  XI ^  for  Z^,  and  the  above  definitions,  one  has, 
from  Equations  (6),  (5),  and  (7): 

N 


/  (Y.  -  Y)  •  XI 

i-  i  '  i 

i-1 

n  71  \X 

;’.xliB  -rf  Ixli 

i-1  \  i  =1 


(8) 


C  -  Y  ■  i-  A 


XI 


i  ’ 


(9) 


i^l 


Equation  (10),  which  represents  the  sum  of  squares  of  Y  residuals,  is 
therefore  equivalent  to  Equation  (2)  but  is  expressed  in  terms  of  only 
the  one  parameter  B. 

If  Q  is  to  be  a  minimum,  the  partial  derivatives  of  Q  with  re- 

k 

spect  to  the  parameters  must  be  zero.  However,  since  A  and  C  are 
determined  from  Equations  (8)  and  (9),  there  is  no  need  to  obtain 
the  partial  derivatives  of  Q  with  respect  to  A  and  C.  From  Equation 
(2 ) ,  one  has : 

N 

^|=2*  L  (Yi  -  A  •  XI®  -  C)  •  (-A  •  XI®  •  In  XI  )  =  0  , 

1=1 

N  N  N 

V  Y.  '  XI®  •  In  XI  -  A  •  '■  XI?  •  In  XI  -  C  •  V  XI®  •  In  XI  =  0. 

_-j  i  r  i  L->  i  l  Lj  i.  i 

1=1  i=l  1*1 

Let  G  represent  the  above  function.  Therefore, 

N  N  N 

G  *  >  Y.  •  XI®  •  In  Xl4  -  A  '  XI?  *  In  XI  -  C  Y  XI®  •  In  XI  . 

2-  i  1  1  <  -  1  1  /  j  1  i 

i-1  i-1  i-1 

(11) 

The  problem  then  becomes  one  of  finding  the  value  of  B  that  makes  G 
zero . 

* 

As  was  stated  previously,  this  is  not  a  sufficient  condition  for 
a  function  that  is  not  1’near  with  respect  to  all  of  its  parameters, 
(bee  Appendix  A. ) 


The  sequence  of  operations  in  -he  computer  program  is  as  follows. 
First,  the  various  suninations  involved  in  Equations  (8),  (9),  md  (11) 
are  obtained  using  B  =  -3.96  (initially).  Then  A  and  C  are  determined 
from  Equations  (8)  and  (9),  respectively.  After  these  calculations 
are  made,  the  value  of  G  is  obtained,  and  its  algebraic  sign  is  noted. 
The  machine  then  steps  the  value  of  B  by  +0.05,  repeats  all  of  the 
summations  and  calculations,  and  checks  the  algebraic  sign  of  G  again. 
This  procedure  is  continued  until  the  algebraic  sign  of  G  is  reversed, 
signifying  that  a  solution  lies  somewhere  between  the  previous  value 
of  B  and  the  value  of  B  at  this  cross-over  point. 

At  this  point,  the  program  begins  an  iterative  operation  in  which 
at  each  cross-over  point  the  incremental  step  is  halved  and  the  direc¬ 
tion  of  advance  is  reversed.  This  iterative  procedure  is  done  as  many 
times  as  desired  to  give  any  degree  of  accuracy  required  for  B.  In 
the  program,  this  procedure  is  repeated  until  the  changes  in  the  abso¬ 
lute  values  of  A,  B,  and  C  from  one  iteration  to  the  next  are  each 

-8 

equal  to,  or  less  than,  10 

The  search  for  roots  continues  to  B  *  -0.01.  After  this  point 
is  reached,  the  program  begins  another  search  starting  at  B  =  +0.01 
and  proceeding  by  increments  of  +0.05  out  to  +3.96,  If  no  solution 
at  all  is  found  within  these  limits,  a  statement  to  this  effect  is 
printed,  and  the  program  continues  on  to  the  next  run.  Any  time  a 
solution  is  found  for  A,  B,  and  C,  the  sum  of  squares  of  Y  residuals 
(Q)  is  determined  and  compared  with  the  corresponding  value  for  the 
previous  solution  (if  there  was  one).  The  solution  that  gives  the 
lowest  sum  of  squares  of  Y  residuals  is  stored  temporarily  for  com¬ 
parison  with  any  future  solution  so  obtained.  In  this  ’/ay,  when  the 
search  is  completed  and  if  there  is  a  solution,  that  solution  will 
generally  represent  the  lowest  sum  oi  squares  of  Y  residuals  -n  the 
region  searched. 

- 

Acknowledgment  is  made  to  Mr.  James  Johnston  (now  at  the  Insti¬ 
tute  for  Defense  Analyses)  for  his  suggestions  in  the  initial  program¬ 
ming  aspects  of  this  problem. 

If  C  is  specified,  then  that  value  is  used  instead  of  calculating 
C  from  Equation  (9). 
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Any  "solution"  found  in  the  specified  range  for  B  represents  a 
solution  for  which  the  partial  derivatives  of  Q  with  vespect  to  the 
parameters  are  zero.  The  Q  value  for  that  solution  is  also  compared 
with  the  Q  values  for  the  end  points  of  B  to  make  sure  that  Q  is  not 
decreasing  to  ~ome  other  minimum  outside  the  range  of  B.  Test  cal¬ 
culation:.  of  ^  as  represented  by  Equation  (10)  indicate  that  Q  in¬ 
creases  smoothly  to  apparent  1;  constant  values  for  large  values  of 
negative  or  positive  B.  As  of  now,  the  author  has  not  been  able  to 
determine  any  requirements  for  unimodality  of  Q  but  has  observed  that 
for  data  applicable  to  this  function,  Q  seems  to  be  unimodal  in  the 
re  .,i'  .  searched.  Even  if  it  is  not,  the  most  minimum  of  those  modes 
will  usually  be  found.  As  stated  before,  one  must  be  aware  that  there 
could  be  other  minima  outside  the  range  searched  which  cannot  be  de- 
’  rmined  by  the  above  method.  However,  this  may  be  relatively  unim¬ 
portant  if  the  "solution"  found  satisfies  the  analyst's  criterion  for 
a  good  fit. 

The  above  limits  on  B  and  the  incremental  stepping  value  of  0.05 
were  chosen  on  the  basis  of  economic  computer  operating  time  and  the 
extent  to  which  the  search  range  should  be  covered  in  order  to  lessen 
the  chances  of  missing  a  root.  Although  two  roots  could  c cr.  ,'ivably 
be  missed  in  the  Incremental  step  of  0.05,  indicating  that  the  G  func¬ 
tion  goes  i.om,  say,  a  positive  to  a  negative  to  a  positive  value  within 
an  interval  of  B  equal  to  0.05,  this  seems  rather  unlikely-  Such  a 
function  would  have  to  behave  extremely  erratically,  and  t<  st  results 
Seem  tc  indicat  2  that  this  function  does  not  gep-rallv  behave  *n  this 
manner . 

Perhaps  it  should  be  noted  that  a  degenerate,  or  trivial,  case 
results  if  B  =  0  or  if  all  Y  values  are  constant  or  all  XI  values  are 
constant.  Any  of  these  conditions  causes  the  numerator  and  denominator 
of  Equation  (8)  to  be  zero  (an  indeterminate  condition).  Of  course, 
this  can  be  seen  from  the  as ymptot ic -power  equation  itself.  Any  of 
the  ’hove  conditions  cause  it  to  reduce  to: 

Y  =  constant . 
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ilBFTC  MIN 

COMMON  A,  Al,  ft,  Bl,  C,  Cl,  0,  f)l  ,  AH,  BH,  CH,  DH,  AN,  AN  l , 

C  AN2 »  OF  1 ,  OF 2 ,  OFT,  LVALUE.  CO,  CV,  R ,  PDEVM,  SEY, 

C  XV,  YV,  YOEVSC,  EA 

COMMON  It,  IDISK,  IERR,  IEV,  JSOLVE,  I  SC l  ,  N,  NFV,  NOTE,  NP 

COMMON  H(4,4),  T I  4  )  ,  S0EV<4),  VMEANI4I,  S(30),  VI201.6), 

C  P0EVI200),  X 1L  (  200 )  ,  X2LI200),  X3LI200*,  YL(200>, 

(  YC1200),  Y0EVI200) 

COMMON  I  SC ( 4  J  ,  JX ( 6 ) ,  FMT(IB),  TITLE!  16),  IND,  IPAGE,  I E  Q- 

C  I ORO ,  NIO,  NRON 

DIMENSION  101(201),  102(201),  Y<201)5  XH201I,  X2(20l), 

C  X 3( 20i ) ,  F  < l ) 

EQUIVALENCE  !FU),  A),  (IOIU),  V(l,l)l,  1102(11*  V<1, 2)). 

C  (  Y  (  1  ) ,  V ( 1 , 3  I  )  ,  iXI(l),  V'(  1,4)),  1X2(1),  V (  1 , 5  )  ) , 

C  1X3(1),  V  (  1,6)  )  ,  (JX<1),  Jl),  (  JX  i  2  ) ,  J2  )  , 

C  (  J  X  (  3  )  ,  J3  )  ,  !jX(4),  J4  )  ,  (  JX  <  5  t  ,  J5  )  ,  (  jX(6)  ,  J6) 

C 

C  MAIN  ROmtinp  nf  • r YE  S '  REGRESSION  ANALYSIS  COMPUTER  PROGRAM 

C 

C  SET  PAGE  NUMBER  TO  t. 

10  IPAGE  «  1 

C  SET  RUN  COUNTER  TO  1. 

NRUN  ■  l 

C  PRINT  MAIN  HEADING  ON  FIRST  PAGt . 

mRITE  (6,  20)  IPAGE 

20  FORMAT  <  1 H 1  / ,  4SX ,  43HCUHVES  REGRESSION  ANALYSIS  COMPUTER  PROGRAM, 
C  27X,  SHPAGE  ,  12  ////  J 

C  STEP  PAGE  NUMBER  B,  l. 

IPAGE  •  IPAGE  ♦  l 

C  U.EAk  common. 

30  00  40  I  >  1,  2706 

HI)  *  0. 

40  CONTINUE 

C  SET  SUBROUTINE  INDICATOR  TO  l. 

I  NO  *  l 
CALL  RE Ai) 

CAi  INPUT 
CALl  check 

r  CHECK  ERROR  DESIGNATOR. 

IF  (  1  ERR  .EU.  I )  GU  TO  1 10 
CALL  SUMS 

go  rn  iso,  60,  jo,  bo,  901,  ieq 

50  CALL  l iNf 
GO  TO  100 
60  CALL  PARA 
SO  TO  100 
70  L  AL  L  PO«R 
GO  TO  100 
80  CALL  ASYM 
GO  TO  100 
90  CALL  EXPO 

100  if  ( IERR  MU.  I )  GO  TO  l 10 
C  At  L  MAT 
CALL  OUT  I 
CALL  OUT  2 

(.  ,T{  P  RUN  COUNTER  BY  1. 
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HO  NRUN  *  NRUN  ♦  i 

RtStT  P AGE  NUMBER  TO  1 
I  PAGE  *  l 
GO  TO  30 

END 


o  o  o  o  o 


H  B  F  T  C  RE  An 


SilSROUT  i  ME  READ 


COMMON 

C 

C 

COMMON 

COMMON 

C 

C 

COMMON 

C 

DIMENSION 

C 

DIMENSION 
EUUI  VALENCE 

C 

C 

C 

DATA 

C 


A,  Al,  8 «  81,  C,  f’,  D,  01,  AH,  BH»  CH,  DH,  AN,  AN  I  , 
AN2,  Dc l ,  0F2,  OFT,  F VALUE ,  CD,  CV,  R,  PDEVM,  S,  , 
XV,  W,  VD  c  V  SQ ,  £  A 

! A ,  10ISK,  SERR,  IFV,  1S0LVE,  I  SC  I  *  N,  NIV,  NOTE,  NR 
H(4,4)f  T(45,  $0EV(4),  VMFANI4),  SI30),  V(20l,6), 

PDF  V(  2001  ,  X 1 L  I  200  )  ,  X2L(2oO),  X3LI200),  YU20G), 

YC  (  200 ) ,  YDEVI200) 

1  SC  I  4 )  f  JX  ( 6  I  ,  FMTI185,  TIUE1I6),  IND,  IPAGE,  IEO. 

I ORD,  NIO,  NRl'M 

101(2011,  ID2W01),  YI20I5,  Xl«20n,  X2(2011, 

X3(201)  ,  Fill 
KX ( 5) ,  K  V ( 6 1 

(Fill,  A),  UOim,  VII,  111,  (102(1),  V(l  ,215, 
i y <  i ) ,  vu,3D,  mm,  v ( i , a 1 1 ,  i x 2 ( i } «  vu.sn, 

(  X  3  (  1 )  ,  VI  1,6)3*  1JX(1),  Jil,  (  JX  (  2  ) ,  J2 )  , 

<  J  X ( 3 ) ,  J3i,  ( JXl 4 1  ,  J4J ,  I JX ( 5  I  ,  J5 ) «  (JXI61,  Jei 
KV(  1),  K Vi ^ )  ,  K V I  3  »  V  K  V  {  4 )  ,  KV(5),  KVC6)  / 

2H ID,  2HY1,  2HX1,  2HX2,  2HX3,  2H  / 


SUBROUTINE  FOR  RFADING  TITLE  CARD,  ORDER  CARD,  AND  VARIABLE-FORMAT 
CARD 


C 


C 


C 

c 

c 

c 


c 


( 


c 


READ  TITLE  CARO. 

10  READ  <5,  20)  I*MT,  IOISK,  IEQ1,  I0RD1,  ISCi,  IA, 

C  (TITLE  II),  I  *  l«  16) 

20  FORMAT  (611,  2X,  16A4) 

CHECK  TITLE  CARD  FOR  TERMINATION  DESIGNATOR  (IFMT  *  °) 

30  IF  ( IFMT  .EC.  9)  CALL  EXIT 
IF  ( I  EG  1  .GF.  IS  iEQ  *  IFQl 
IF  U0RD1  .GT.  0)  I  ORD  -  1  -  (KRD1/2) 

PRINT  TITLE  AND  PAGE  NUMBER  ON  FIRST  PAGE 
40  if  (NPUN  .G  ' .  1)  GO  TO  55 

MRIIE  (6,  50)  (TITLE!!),  I  *  1,  16) 

50  FORMAT  ( 1H  /  ICX,  16A4  ///  ) 

GO  TO  65 

55  WRITE  (6,  60)  ( T ! Tl  E ( I ) ,  I  *  l,  16),  IPAGE 
60  FORMAT  (  1H 1  /  10X,  16A4,  41X,  5HPAGE  ,  12  ///  ) 

STEP  PAGE  NUMBER  BY  I. 

IPAGE  «  IPAGE  *  1 
PRINT  SUBHEA01NG. 

64  CALL  PRINT 

TEST  WHETHER  TITLE  CARD  FOR  FIRST  RUN  CONTAINS  THE  V AR I A8L E-EORM AT 
INDICATOR  IN  COLUMN  I. 

IF  ( NRUN  .GT.  1)  GO  TO  90 
IF  <  IFMT  .GT.  0)  GO  TO  IOC 
WRITE  ERROR  MESSAGE. 

TO  WRI1F  (6,  80) 

80  FOR 4  A  T  UHO//,  IQX,  U8HTHE  VARIABLE-FORMAT  INDICATOR  HAS  NOT  BEEN 
(.  ENTERED  IN  COLUMN  l  OF  THE  FIRST  TITLE  CARD.  THIS  JOB  HAS  BEEN  T 
TERMINATED. ) 

CALL  EXIT 

TEST  WHETHER  ORDER  CARD  AND  FORMAT  CARD  ARE  TO  BE  READ  FOR  THIS 
RUN. 

90  IF  (  IFMT  .EG.  0)  GO  TO  170 


C.  RE  An  DROFf’  CARD. 

100  READ  { 5 1  KO)  UXU).  l  *  it  5>,  NIDi 
110  FORMAT  (5 »A2t  1X}«  1 1 / 

IF  (NIDI  .GT.  01  NiO  -  MIDI 
I  «  \ 

JMAX  »  1 
00  ISO  K  *  1,  5 
DO  120  J  -  l,  6 

IF  (KXOU  .ED.  KVjj))  SO  Tf  30 
120  CONTINUE 

C  WRITE  ERROR  MESSAGE. 

WAITE  I  6 «  125  J 

125  FORMAT  < 1H0//  10X,  98HTHERE  IS  AW  ERROR  IN  THE  VARIABLE-FORMAT 
C  INDEX  (SECONO  INPUT  CARO!.  THIS  JOB  WAS  BEEN  TERMINATED.  ! 

CALL  EXIT 

no  if  (j  .ne.  n  go  ro  no 
jxu )  -  i 

JXUU)  >  2 
I  «  S  ♦  2 
GO  TO  145 

140  jxm  «  j  ♦  i 

I  *  I  ♦  1 

145  IF  (J  .EQ.  6)  GO  TO  150 
IF  UMAX  .IT.  Jl  JMAX  -  J 
150  CONTINUE 
C  READ  FORMAT  CARD 

READ  (5,  It  >1  (FMTI II,  I  •  1,  18) 

160  FORMAT  (1BA4) 

C  NUMBER  OF  INDEPENDENT  VARIABLES 

ITO  MI V  «  JMAX  -  2 

C  NUMBER  OF  PARAMETERS  TO  BE  SOLVED  (NP) 

NP  •  M!V  *  l  -  I A 

IF  UEQ  «EQ.  2  .OR.  IEO  .EQ.  41  NP  «  N'»  *■  1 
IF  1 1 A  .CO.  0  . AND.  ISC  I  .ED.  0)  RETURN 
C  READ  SCALE  FACTORS  AND/OR  VALUE  OF  V  INTERCEPT. 

IF  I  IEO  .NC.  4)  RE  AO  (5.  ISO)  IlSCdi,  I  -  1,  41,  A 
IF  (IEQ  .EQ.  4)  READ  15,  160)  (1SC(U»  I  -  1,  4),  C 
160  FORMAT  (412,  tR,  F10.0) 

RETURN 

ENO 


t'HFTC  P3  |M 


C 

c 

c 


SUBROUTINE  PPiNf 

COP!?  ON  A  ,  A  1  ,  8 ,  i3l,  Cc 

ANA  nf  L  ,  OF?,  I'M 


,  D,  01, 

►  Vi*  s ifc  , 


AH,  PH,  '  H,  DH,  AN,  AN  1 , 
CO,  Cw;  R «  ?»ncVWs  $fvt 


C 


COMMON 

COMMON 

C 

C 

COMMON 

C 

DIMENSION 

L 

tOUl VALENCE 

C 

c 


W.V,  vv,  YDCVSu,  E A 

IA,  1D1SX,  IERR,  I  F  V ,  ISOLVE,  ISClt  N,  NlV,  NOTE,  NP 
H(4,4),  T  (  4  )  ,  SDEVI4),  VMEAN(4),  S(3G»,  V»20i,6), 

POE V ( 2001 ,  X1U200),  X2LI200),  X3LI200),  YLI20C*. 

YC ( 200 ) ,  YOE V ( 200) 

l  SC ( 4 }  »  JX  (  6 )  ,  F  M  T  !  i  R  J  ,  IULtllM,  INO,  I  PAGE ,  IEO, 
10RO,  NIO,  NRUN 

(01*201),  102(201),  Y  (  2  0  1 )  ,  Xl(.f»l)»  X2(?0!S 
X  3( 20  l  J  ,  F  <  1) 

( f ( i ) ,  a),  noun,  vu, m,  iidoiii,  vn,2n, 

(Yd),  VII,))),  IXl(l),  V  (  1 , 4  )  )  ,  (  X2(  1 )  ,  V  '  1*5)  )  , 
1X3111*  vll,6!i,  Ullll,  Jl),  (  JX  (  2 )  *  J2), 

( JX ( 35  «  J3 ) ,  ( JX ( 4 ) ,  J 4 ) ,  *JX<5>,  J5),  (JX(6>,  J6) 


SU8P0UT 1 NE  FOR  PRINTING  SUBHEADINGS 


10  GO  TO  (20,  40,  60,  80,  100),  I E 0 

20  WRITE  (6,  30) 

30  FORMAT  (LHC/,  3faX,  59HL  i  NEAR  REGRESSION  -  T  *  A  ♦  IB  Ml)  HC  ♦ 
C  X2)  ♦  I D  *  X 3)  /  ) 

RETURN 

40  WRITE  (6,  50) 

50  FORMAT  (  1H0/ ,  39X,  5 2HPAR ABQl  I C  REGRESSION  —  '/  *  U  ID  *  XI  ♦ 

C  (C  «  X**2)  /  ) 

RETURN 

60  WRITE  (6,  TO) 

70  FORMAT  (  1H0/  37X,  55HPOWER  REGRESSION  —  Y  =  A  *  (Xl**8?  *  (X2<'AC) 

*  #  (X3»*D)  /  ) 

RE  TUN 

an  WRITE  (6,  90) 

90  FORMAT  (  l HO/ ,  39X ,  S lHASYMP TO? IC-POWER  REGRESSION  —  Y  *  (A  * 

<  (X*ABH  4  C  /  ) 

RETURN 

IOC  WRITE  (6,  110) 

IIC  FORMAT  ( 1  HO / ,  42X,  46HE XPONF NT  I AL  REGRESSION  —  Y  *  EXPIA  ♦ 

♦  ( t)  *  X  H  •  • 

«f  TURN 


(  NO 


l 

i 

t 

i 

i 

j 

4 

1 

| 


5 


1 


I 


o  r.  o 


1 1 


l  NPUT 


jUBP  DU  T ! Nfc  INPUT 


OMMCN 


L 

C 

COMMON 

COMMON 

r 

■u 

c 

COMMON 

C 


DIMENSION 

C 

DIMENSION 

equivalence 

C 

C 

C 

DATA 


A.  A l f  fc ,  R '  ,  Ct  Cl,  0,  01,  AH,  BH.  CH,  DH,  AN,  491, 
AN? ,  OEl,  0F2,  OF  T ,  F VALUE •  CD,  CV,  R,  POEVM,  Sfv, 
XV,  YV,  YOEVSQ,  FA 

IAf  '.DISK,  I  ERR ,  ifV,  ISOLVE,  I  SC  l  •  N,  NIV,  NOTt,  NP 
H I  4  »  4 ) ,  T ( 4 ) ,  SDEV(4),  VME  AN  I  4  I ,  SI30),  V(20l,6>, 

PDF  V  { 200 )  ,  X  IL  (  200 )  ,  X2LI200),  X3LI200),  YL1200). 

VC  I  2001  ,  YDEVI 200) 

I  SC ( 4) ,  JX ( 6 ) ,  FMT 118),  TIiLEII6)«  IND,  IPAGE,  IEQ, 

I ORD,  N I  «)»  NRUN 

(01(201),  102(201),  v i 201 ) »  X 1 ( 201 ) ,  X2(20l), 

X  3 ( 20 1 ) ,  F ( 1 ) 

V0ATAI20) 

(Fit),  a ) ,  (IDicn,  Vll.li),  UD2(l),  V(l,2)), 

( Y ( l ) ,  V( 1 , 3 ) ) i  (  X  1 1  1 )  *  VI  1,4)1,  ( X2I  1 ) *  VII, 5)), 
1X3(1),  VII,  6)),  (Jill),  JU,  (  J  X I  2 )  ,  J2)  * 

IJXI3),  J3 ) ,  ( JX ( 4 ) ,  J4),  IJX(5),  Jal,  tJ2U6>,  J6) 
BLANK,  (BLANK/4H  ,  4H  / 


SUBROUTINE  FOR  READING  IN  DATA  ANO  WRITING  DATA  ONTO  UTILITY  DISK 


10  IF  I  I D I SK  .EQ.  0)  GO  TO  60 
REWIND  4 

C  SET  INPUT  TAPE  NUMBER 

IM  *  4 

IF  « (DISK  .EQ.  2)  GO  TO  70 
00  SO  I  *  1,  201 

C  READ  INPUT  DATA  AS  ALPHANUMERIC  DATA 

READ  (5,  20)  (VOATAIJ),  J  *  l,  201 
20  FORMAT  (20A4) 

C  WRITE  INPUT  DATA  ONTO  UTILITY  DISK. 

WRITE  (4,  20)  IVDATA(J),  J  «  I,  20) 

C  CHECK  FOR  BLANK  CARO. 

00  30  K  -  l,  20 

IF  (VOATA(K)  .NE.  BLANK)  GO  TO  50 
30  CONTINUE 
REWINO  4 
40  GO  TO  70 
50  CONTINUE 
KEWINO  4 
GO  TO  70 

C  SET  INPUT  TAPE  NUMBER. 

60  IM  *  5 

C  READ  INPUT  DATA  FROM  EITHER  CAROS  (IDISK  =  0)  OR  FROM  UTILITY 

C  OISK  (IOISK  -  1). 

70  (READ  *  0 

DO  80  I  *  1,  6 

If  (JX!  I  >  .EO.  I )  I  ‘AO  *  l 
80  CONTINUE 

C  SET  IDENTIFIERS  TO  BLANK 

DO  BS  I  =  l,  200 
IC1 (I  )  *  IBIANK 
102(1)  *  I  BLANK 
85  CONTINUE 

IF  (  I" *  AO  .tV.  I )  GO  «0  90 


GO  Tn  t  too,  1  10,  1601 ,  *  IV 

90  ! f  t n i o  .eg.  n  go  to  28o 

GO  TO  ( 190,  220,  2505 ,  Nl V 
100  00  120  !  *  1,  201 
110  REAO  (IN,  FNT )  VII, Jl),  VI  I  *  J2  l 

IF  (A8SIYIU)  ♦  ABSIXl(II)  .£Q.  0.)  GO  to  580 
IF  (Yin  .  EQ.  0.  .OR.  Kill)  -£Q.  0.)  GO  TO  110 
120  CONTINUE 
GO  TO  590 

130  00  150  I  *  1 ,  201 

150  READ  UN,  FHD  VII, Jll,  VII, J2J,  VII, J3) 

IF  lABSlYlDI  *  ABSIXIID)  >  ABS(X2im  -EQ.  0.)  GO  TO  580 
IF  (YU)  .EQ.  0.  .OR.  XIII)  .EQ.  0.  .OR.  X2I1)  .EQ.  0.)  GO  TO  140 
150  CONTINUE 
CO  TO  590 

160  DO  180  1  «  1,  201 

170  READ  (IN,  FNT )  VII, Jl),  VII, J2),  VII, J3),  VII, J4) 

IF  (ABSIYim  ♦  ABSIXIID)  ♦  ABSIX2U))  ♦  A6SIX3U))  .EQ.  0.) 

C  GO  TO  580 

IF  |  Y  (  I  )  .EQ.  0.  .OR.  XIII)  .EQ.  0.  .OR.  X2U)  .EQ.  0.  .OR. 

C  X3I  I)  • F O  j.)  GO  TO  170 
180  CONTINUE 
GO  TO  590 

190  00  210  I  *  1.  201 

200  READ  (IN,  FNT)  VII, Jl),  VII, J2),  VII, J3),  VII, J4) 

IF  (ABSIYim  ♦  ABSIXIID)  .EQ.  0.)  GO  TO  580 
IF  (YID  .EQ.  0.  .OR.  XIII)  .EQ.  0.)  GO  TO  200 
2  1C  CONTINUE 
GO  TO  590 

220  UO  240  I  *  1,  201 

2 30  READ  (IN,  FNT)  VII, Jl),  VII, J2),  VII, J3),  VII, J4) 

IF  (ABSIYID)  ♦  ARSIXl(D)  ♦  ABS I X2 1 1)  I  .EQ.  0.1  GO  TO  580 
If  (YID  . EQ.  0.  .OR.  XIII)  .EQ.  0.  .OR.  X2I1I  .EQ.  0.)  GO  TO  230 

2iO  CONTINUE 
GO  TO  590 

250  DO  270  I  «  l,  201 

260  RE  AO  UN,  FNT)  VU,Jl),  VII, J2),  VII,  J3>,  VII. J4),  V 1 1  ,  J5 ) ,  VI  I  ,  J6 
IF  (ABSIYID)  ♦  ABSIXIID)  ♦  ABS  I X2 1  I  ) )  ♦  A8SIX3U))  .EG.  0.) 

C  GO  TO  580 

If  IYIII  .EQ.  0.  .OR.  XIII)  .EQ.  0.  .OR.  X2U)  .EQ.  0.  .OR. 
l  X  3(1 )  .EQ.  0. )  GO  TO  260 
270  CONTINUE 
280  00  290  X  *  1,  5 

IE  ( JX ( K )  .EQ.  1)  GO  TO  300 
290  CONTINUE 
GO  TO  610 

300  GO  TO  (310,  380,  470),  NlV 
HO  GO  TO  I  320,  340,  360)  ,  K 
120  00  330  I  «  1,  201 

325  Rf  AO  UN,  FNT)  VUtJU.  VU,J3),  VII, J4) 

IF  (ABSIYim  ♦  ABSIXIID)  .EQ.  0.)  GO  TO  500 
IF  (YU)  .EQ.  0.  .OR.  XIII)  .Es..  0.)  GO  TO  325 
330  CONTINUE 

GO  TO  590 

140  00  350  I  -  1,  201 

145  READ  UN,  FNT)  VII,  Jl),  V(I,J2),  VII, J4) 


If  (ABSUim  *  «^!Uim  .to.  0.)  go  to  5«o 
(F  } r( I J  .EQ.  0.  .OR.  X1(I)  .EQ.  0.)  GO  TO  345 
350  CONTINUE 
GO  TO  590 

3<S0  00  3  70  I  *  i,  ?0i 

365  RE  AO  UN,  FNT)  Vd.JU,  Vi  f  •  J2  >  •  V  (  E  »  J  3  3 

IF  (ABS(YU))  ♦  ABSUlim  .EQ.  O.J  GO  TO  580 
IE  1 Y 1  I )  .EO.  0.  .OR.  XIII)  . EO.  0. )  GO  TO  365 
370  CONI INUE 
GO  TO  590 

380  GO  TO  (390,  410,  430,  450),  K 
390  00  400  I  *  1,  201 

395  READ  UN,  FNT)  VIUJ1),  VII. J3i,  VII, J4),  VII,  J5> 

IF  (ABSmiH  ♦  ABSUlim  ♦  ABS(X2U))  .EO.  0.)  GO  TO  *80 
IF  lYlU  .EO.  0.  .OR.  Klin  .EQ.  0.  .OR.  X2I  I)  .EO.  0.  >  TO  395 
400  CONTINUE 
GO  TO  590 

410  00  420  I  «  1,  201 

415  READ  (IN,  FNT)  V(I»J1I,  Vii,J2),  VU.J4),  VII, J5) 

IF  (ABSIYII))  ♦  ASS (XIII)}  ♦  AB5I  X2U  )  )  .EQ.  0.)  GO  TO  5B0 
IF  (YU)  .EQ.  0.  .OR.  XIII)  .EQ.  0.  .OR.  X2UJ  .EQ.  0.)  GO  TO  K5 
420  CONTINUE 
GO  TO  590 

430  DO  440  I  »  l,  201 

435  READ  (IN,  FNT  I  VU,J1),  VU.J2I,  V(I,J3),  VII, J5» 

IF  (ABSIY III)  ♦  ABSUlim  ♦  ABSU2I1II  .EQ.  0.)  GO  TO  5B0 
IF  ( ri I )  .EQ.  0.  .OR.  XIII)  .EQ.  0.  .OR.  X2II)  .EQ.  0.)  GO  TC  *35 
440  CONTINUE 
GO  TO  590 

450  CO  UO  i  -  i,  cui 

455  READ  UN.  FMT)  VU,J1>,  VU.J2),  VII, J3),  VII, J4V 

IF  (ABSmm  ♦  ABSUlim  ♦  AB$(X2im  .EQ.  0.)  GO  TO  580 
IF  tYUl  .EQ.  0.  .OR.  kill!  .EQ.  0.  .OR.  X2II)  .EQ.  O.I  GO  TO  455 
460  CONTINUE 
GO  TO  590 

470  Gu  TO  1480,  500,  520,  540,  560),  X 
480  00  490  I  •  I,  201 

485  READ  <IN,  FNT  t  VII,  JO,  VU.J3I,  VII, J4),  V  U  ,  J  5  )  .  VII,  J6> 

IF  (ABSIYU)I  ♦  ABSUIU))  ♦  ABSU2III)  ♦  ABSUXII)  .EQ.  0.1 

C  GO  TO  580 

IF  (YII)  .FO.  0.  .OR.  XKI)  .EQ.  0.  .OR.  X2U)  .EQ.  0.  .OR. 

C  X3U)  .EQ.  0.  )  GO  TO  485 
490  CONTINUE 
GO  TO  590 

500  DO  510  I  •  1,  201 

505  READ  (1N,  FNT)  V(I,Jl),  V(l,J2),  VII, J4),  VII, J5 > ,  VII, J6) 

IF  U6SIYU))  ♦  ABSUlim  *  ABSU2ID)  ♦  ABSU31II)  .EQ.  0.) 

C  GO  TO  580 

IF  (YU)  .EQ.  0.  .OR.  XIII)  .EQ..  0.  .OR.  X2  (  I  I  .EQ.  0.  .OR. 

C  XKI)  .EQ.  0.  )  GO  TO  505 
510  CONTINUE 
GO  TO  590 

520  00  530  I  «  I,  201 

525  READ  (IN,  FNT)  VII, Jl),  VII. J 2),  V(I,J3»,  VII, J5I,  V(I,J6l 

IF  (ABS(YU))  ♦  ABSUlim  ♦  ABS(X2U))  ABSUHIU  .EQ.  0.) 

C  GO  TO  580 


1 f  <MI1  .  u  .  o . 


x  3 (  .  )  0. 


.  t  r 


o .  )  (,  ( 


<  i  (  I  i  .  f 

s/s 


;  \a 


*  /  (  !  )  .  f  Q 


.  i',9 


53C  L ON  F  I  Nut 

on  TO  S90 

‘jtn  nn  kfcn  i  s 

“II  ,  J?)  ,  VII. JO,  VII, JO.  VII,  J6I 

c  00  *0  Uo  *  *"SUI""  *  •  MSIOini  .to.  0.1 

c 1  * » ! 1 5 ' !to!°o.i'Gc0[o  w"  -eu-  °-  -nR- ,?m  °- 

5  5C  CONI  I  NUE 
or  TO  5 90 

560  00  570  I  *  1,  201 

565  READ  UM,  FMTI  Vd.Jl),  VU.J2),  VII. JiJ.  V(I.J4>,  V(l  J5J 

C  "niello1"  4  *BS,xmn  ♦  ♦  ABS(x3<m  !eq.  o. , 

<.,x-.!n'!tS!‘i;.?-on°;o  -E0-  °-  -0R- -te-  »•  ■0B- 

5  70  CONTINUE 
GO  TO  590 

5«0  IF  (I0ISX  ,NE.  0 )  REWIND  4 

SET  N  EOUAL  TO  NUMBER  OF  0A'">  pr,f*js 
N  *  I  -  1  * 

FLOAT  N. 

AN  *=  N 

NUMBER  OF  OATA  POINTS  LESS  » 

AN  1  *  AN  -  1.0 

number  of  DATA  POINTS  LESS  2 
AN?  *  AN  -  ?.0 
B  F  T  |  |U  »■ 

5  90  WK|  TE  (6,  6l-0) 

,,a"i  if"''  <">MNU»et«  St  INPUT  0«T*  POINTS  »»s  f.cttoto 

t«um  ''n<”‘  Tf"S  J,,p  M*s  8">  Tf««IMTE0.l 

MO  WRITE  (6,  620 ) 

620  FORMAT  (  ,HG//  IQX,  111HTHEKE  IS  AN  ERROR  ASSOCIATED  v I T  N  I Hf 

CmIn1teM)UM"K1RM4T  i,fcSlGNM0R  ‘0KPf«  l-ARO).  THIS  JOB  HAS  BEEN  TER 

'.All  EXIT 
t  NO 


!  8  E  r  ,  C  Ht  L  f 

';!,ProHM  .  hKk 

A .  At,  rt,  H  t  ,  C,  Li,  0,  Dl.  am,  Rh,  Cm,  pm,  an,  AN  l  * 
C  AN?,  Uf  I  ,  C  F  2  »  CFf,  ivnut  ,  CD.  CV.  R.  P(>FVM,  SfY, 

(,  t  V  ,  Y  V  ,  Y  D  E  V  S  w  .  E  A 

COMMON  (A,  [DISK,  1FRH,  IfV,  ISOlVf,  (  SC  i  .  N ,  N[v,  NOTE-  NP 

COMMON  Ht4,4),  T  (  <,  :  ,  SDEV(M,  VME  AN ( A )  ,  S<30),  V!20l,fc), 

L  P0EVI200),  x 1 1  ( ?  00  )  ,  X  2L ( 200  ? ,  X3H200),  vl(200), 

C  YC(?CO),  Y  0  E  V  f  2  0  0 ) 

COMMON  I  SC  (  A  )  ,  J  X  I  6 )  ,  f  M  T  t  i  8  )  ,  MUEI16),  l  NO  ,  (  PAGE  ,  lEO, 

C  I  CRD  *  N10,  NRUN 

DIMENSION  101(201)?  102(201),  v ( 20 1 )  ,  Xl(?Ol),  X2(?01), 

C  X3( 2C1) .  F (!) 

DIMENSION  R  5C ( A ) 

EQUIVALENCE  ( F ( 1 ) ,  A),  (101(1),  V (  1,1)),  (  1 D? i  i  ) »  V (  1*2)), 

r.  ( r  <  i ) ,  v  1 1 , 3  n  ,  ( x  i  m  ,  v  ( 1 ,  a  n  ,  <  x  2  ( 1 ) ,  v  1 1 ,  s  >  j , 

C  1X3(11,  v(  1 , 6 )  ) ,  (JX(l),  Jll,  ( J  X ( 2  )  *  J  2  )  « 

L  (  J  X  !  3  )  ,  J)),  (  J  X  (  A  )  ,  JA),  (JX(S),  J  *>  I  ,  (  J  X  (  6  )  ,  J  6  ) 

SUBROUTINE  FUR  CHECMNG  0  A I  A  FCR  ERRORS  AND  SCALING  DATA 

CALCULATE  SCALE  FACTORS. 

10  DO  2C  I  *  1,  A 
HSC  f  I )  =  1.0 

IF  (iscm  .  NE„  0)  RSCU)  *  I  0  -  0  *  »  (  I  SC  (  I  )  J 
2 C  CONTINUE 

3C  CO  TO  (AO,  60,  HOI,  MV 

AO  CO  SC  I  *  t,  N 

If  (  Y  (  J  )  .IF.  0.  .OR.  X 1  (  I  )  .LF.  C.  )  '0  TO  170 

DO  AS  j  *  3,  6 

v(!,J)  *  VII,J)  *  RSCU -2) 

*.5  CONTINUE 
SO  CONTINUE 
GO  TO  120 


6  0 

DO  70  ! 

*  1 . 

N 

if  (Yin 

.LE 

,  0. 

.or.  x  n  n 

f  .  ■'  . 

; ;  3 

* 

X  ?  :  I 

(  .  L  E  . 

DO  6  5  J 

-  3, 

6 

V( I , J)  * 

V(  1 

,  J  1 

♦  RSlI.I-2) 

6S 

CONTINUE 

7  C 

CONT INUt 
GO  TC  l« 

V.' 

8  C 

Lfl  9C  I 

1  l , 

N 

IF  (  Y(  1  ) 

.IF 

.  C. 

.or.  XII!) 

.A.  C.  . 

* 

x  ? !  ! 

>  .  L  E  . 

C  X  3  (  I  )  . 

IE. 

0.  ) 

GO  ti  I/O 

CO  8  S  J 

■  3 , 

6 

vino  ■ 

v  i  1 

,  J  ) 

*  R  SC  (  «-2  ) 

8S 

CONT [Nut 

RC 

CCN  T  I  N L  (' 

120 

I  f  !  I  SC  ( 

if  • 

Nf  . 

0)  MR! TE  ( 6  . 

1  30  i  K SC 

i  \ 

! 

If  (  1  SC  (  2  5  . 

NE  . 

0  )  m  R  I  T  E  (  6  , 

I  A 0  i  R  Sl 

{  <L 

) 

IF  (  1  SC ( 

> .  , 

Nf  , 

C)  MHlTf  (ft, 

ISO)  8  SC 

(  3  ) 

1 F  !  ISO  i 

A  )  . 

Nf  . 

0  )  MR  ?  Tf  |  fjj 

Ifti'l  RSl 

1  * 

} 

S  5  0 

f  OR*  a  T  ( 

1  H  , 

OX  , 

2 :  v  SCALE 

E  AC  T  i  a 

,  18. 

1  ) 

•  *,  i 

*  !  R  “  A  T  { 

1  M  , 

R*  , 

2 OH  X 1  S C  A i  E 

U.  ;  OR 

,  f  H  . 

1  ) 

1  s 

FORMAT  ! 

1  M  , 

8  *  , 

2 1  I  •  <  2  SCALE 

f  M  If,: 

,  Eh. 

l  ) 

,  ' 

s  R  “  A  i  ; 

1  M  , 

NX  , 

7(  MX  '  SCAtE 

FAC  I  0  R 

%  f  M  • 

l  ) 

~  s  K  - 


Rt TURN 

C  ERROR  MESSAGE 

1  TO  WRITE  (6,  180  I 

ISC  FORMAT  (  1  HO / /  10X,  81HA  ZERO  OR  NEGATIVE  VALUE  EXISTS  IN  THE  INPUT 
C  OATA  ( E C°  EXAMPLE,  OATA  CARO  NUMBER  ,  IT.33H5,  THIS  RUN  HAS  BEtN 

G  TERMINATED.  ) 

I  ERR  *  1 
RE  TURN 
END 


II8FTC 


SUMS 

SUBROUT l NE 

COMMON 


COMMON 

COMMON 


COMMON 


dimension 


EQUIVALENCE 


Sums 

A,  Al,  B,  Bis  C,  Cl,  0,  01,  AH,  BH,  OH,  OH,  AN,  AN l 
AN  2,  DEL,  CF  2 ,  OFT,  FV  ALUf  ,  CD,  CV,  R,  PDEVM,  St  Y  , 
XV,  YV,  YDEVSQ,  t  A 

IA,  IDISK,  !ERH,  I F  V ,  ISOLVE,  I  SC  l  ,  N •  NIV,  NOTE, 

H  I  4  ,  A  1  ,  T  (  4  )  *  SD6  V  (  4  5  ,  VME  AN  I  4  )  ,  SHO),  VI  201 ,6  »  , 

PDEV(200),  XILI20U),  X2i.!200),  X  III  200 )  ,  YLi200), 
YC12C0),  Yrev<200) 

ISCI4),  J  X ( 6 ) ,  F  *  T {  10)  ,  TITLE! 16)  -  IND,  (PAGE,  IE 
10RP.  NID,  NPUN 

i  C  i  i  2  0 1  )  «  102(201 ) ,  r  <  20  1  )  ,  <1(201),  *2(201), 

X3( 201 ) ,  F ( 1 1 

|F(1(,  .*),  I  IDUt),  VI  1,11),  (102(1),  V<  1,2)1, 
(Y(t),  V(l,3»),  (XUl),  V  (  1 , 4 )  )  ,  (X2(i),  V  I  1  »  5  )  )  « 

( X  3 (  1  )  ,  V I  1 , M  1  ,  ( J  *  ;  1  )  ,  J  1)  ,  ! J X ( 2  )  ,  J  2  )  , 

I  JX  (  3  )  ,  JO,  (  J  X  !  *»  )  ,  J  4  )  ,  (JXIS5,  JS),  (  J  X  (  6  )  ,  Jfc 


I  2  AGE  , 


V(  1,2) ) , 
.  )  .  V<  1,5 
J?  )  , 

,  ( JX ( 6)  , 


SUHROUT l Nt  FOR  ORDERING  DATA  AND  OBTAINING  VARIOUS  SUMS  AND 
STANDARD  DEVIATIONS  CF  INPUT  DATA 


IF  (  I  OR  t  .EQ.  0)  GO  TO  «,0 

ORDER  The  DATA  FROM  EC*  TO  HIGH  VALUES 

NK  *  N  ~  ! 

DO  40  I  =  l,  NK 

(N  -  I  ♦  l 

DO  30  J  *  IN,  N 

IF  *  Y (  1  )  ,LE.  Y ( J  )  )  GO  TC  30 

CO  20  K  »  t,  6 

TEMP  •  VII.XI 

VI  i  ,K)  «  V! J,K) 

VIJ.M  -  TEMP 
CON»  INUE 
CONTINUE 
CONTINUE 


08 T A  IN 

VAR IOUS 

sums  c 

F 

INPUT  D 

oo  r 

I  *  1,  N 

YL  (  I) 

•  ALOGI Y I : ) ) 

xil  (  n 

-  ALOGI XI  I  1)) 

Also  • 

XIII)  • 

XIII) 

ST  l  i  - 

sin  ♦ 

Y  I  I  ) 

S  (  2  )  • 

S  I  2  )  ♦ 

X  l  I  1  ) 

S(  3)  • 

SI  31  ♦ 

X2I  1  ) 

S  (  4  1  - 

S  I  4  )  ♦ 

X  3  I  I  ) 

SIS)  • 

SIS)  ♦ 

XI  sc 

S  I  6  )  * 

S  I  6  )  ♦ 

I XI  m 

• 

x  2  ( 1 1 ; 

SI?)  - 

SIM  ♦ 

mil) 

0 

X  3  I  I  )  ) 

SI  8  )  • 

S<  8  )  ♦ 

i  x  2 1 ;  i 

• 

X  2  1  l  )  ) 

S(<*  )  * 

SIN!  ♦ 

1X2(1  ! 

0 

X  3  1  l  )  ) 

S(  10) 

»  SI  10) 

♦  (X3II 

) 

*  X  J  (  !  ) 

S(  111 

-  Sill) 

♦  mm 

♦  V  I  I  »  ) 

S(  12) 

«  S  1  1  2  ) 

♦  1  X .'  (  I 

) 

*  Y  !  I  )  ) 

S(  l  3) 

-  Mill 

♦  (  *  3  (  1 

) 

♦  Y  (  |  )  ) 

IF  (  IEC  .EQ.  1  )  GO  rn 

6  0 

SI  14) 

■  SI  14  ) 

♦  Yl  I  1  ) 

SI  IS) 

•  sim 

♦  X  LL  !  i 

) 

Si  16)  -  $(  If. )  *  {*!(!)  ♦  YLii)] 

s !  17}  =  s<i7»  (mm  ♦  xum) 

IF  (  IEO  ,NE .  2)  GC  TC  6 C 

sub)  -  si  >8)  ♦  c x i s g  *  xmn 

si  m  *  si  19)  4  ms**  *  *isn 

Si  20 )  =  Si  20 i  ♦  ( XI  SO  *  Y  (  i  )  ) 

6C  CONTINUE 

CALC Ut  A  TF  MEANS  CF  INPUT  DATA. 
i;u  70!  -  l  *  A 
V  M  E  A  N I  n  *  Sin  /AN 

70  CONTINUE 

CALCULATE  sums  ff  The  i  N  i;i  i  u  A  T  A  ABCu  T  ThEIR  MEANS. 
CD  90  K  ^  i  ,  A 
DO  80  i  =  U  N 

<(*♦20)  *  S{K*2r')  *  i(v<i,K42J  ~  VMEAN(K))**?) 

HP  CONTINUE 
9C  CONTINUE 

CALCULATE  ST  AND ARC  DEVIATIONS  OF  INPUT  DAT  A„ 

DO  ICC  K  *  l, 

SDfcViK)  -  SORT ( S(k  >20}/ANl ) 

TOC  CONTINUE 
RETURN 
END 


oooooooo 
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$  !  8  F  T  C  LINE 


SUBROUTINE 

COMMON 

r 

C 

COMMON 

COMMON 


C 


c 

c 

C 

C 

c 

c 


COMMON 

DIMENSION 

EQUIVALENCE 


LINE 

*r  Al,  8 1  81,  C,  CL,  D,  Dl,  AH,  8H ,  CH,  OH,  AN,  AN  1  , 
AN2 ,  CF1,  DF2,  OFT,  FVALUE,  CD,  CV,  R,  PDEVM,  SET, 
XV,  VV,  YDfVSg,  EA 

IA,  IDISK.  I  £  R  G  •  I F  V ,  ISOLVE,  I  SC  1  ,  N,  NIV,  NOTE,  NP 
H(A,A),  T  (  '  )  »  SDEV  I  A '  ,  YME  AN ( A  » ,  S < 30 >  ,  Vi20i»6>, 
PDEV { 200) t  X IL ( 200 ) ,  X2L(2GO),  X3LI200),  YU2005, 
YCI2CCJ,  Y0EVI200) 

ISC?A),  JX  <  6 ) ,  FMTiio),  TJTLEII6),  iND,  {PAGE,  IEQ, 

I ORD ,  NID,  NRUN 

101(2011*  102(201),  y ( 20 1 ) ,  X 1  (  201 ) ,  X2I201), 
X3I201),  F  (  l } 

(Ell),  A),  (  I  D  1  (  i )  ,  VU.D),  I  I  D2  (  1 ) ,  V  (1 »  2 ) ) , 

(YU),  VU,3>>,  SXKU,  V  t  i  ,  A  )  )  ,  IX2U).  Vii,5))( 
IX3I1I,  V !  1,6)},  (JXIU,  JU,  (JXI 2),  J2)  , 

(  JX  {  3  )  ,  J  3  )  ,  (JX(4».  i*{,  (  J  X  (  5  )  ,  J5)  ,  jot 


SU1R0UT I N£  FOR  DETERMINING  LEAST-SQUARE  S  SOLUTIONS  OF  PARAMETERS 
FOR  LINEAR  EQUATIONS  OF  FORM 

Y  *  k  ♦  (B*xl)  ♦  ( C  *X  2  )  ♦  ( D*X  3 ) , 


WHERE  A  MAY  BE  SPECIFIED 


C  DETERMINE  VALUES  FOR  ^EFFICIENTS  OF  LINEAR  EQUATIONS  TO  BE 

C  SOLVED. 


1C  IF  (IA 

•  i 

EQ.  1) 

h<  l,  n 

at 

AN 

HI  1,2) 

X 

SI2) 

HU, 33 

= 

S<  35 

HU, A) 

* 

S  1  A  1 

H  (  2 , 2 1 

a 

SI  5 ) 

H  (  2  »  3 ) 

X 

S  ( 6  5 

Hi  2 , A) 

at 

S(  7) 

HI  3,33 

a 

S  ( 8 ) 

HI  3, A) 

= 

$(  9) 

H  (  A  ,  A  ) 

* 

S(  10) 

H(2,  1) 

x 

HU,  2) 

HI  3,  l) 

9 

HU,  3) 

H( 3,2) 

X 

H  (  2 »  3 ) 

H< A,  l) 

X 

HU, A) 

HI  A, 2) 

X 

H  (  2  ,  A ) 

HI  A, 3) 

3 

HI  3, A) 

1111  « 

111) 

T  (  2  )  » 

sun 

T  I  3  i  * 

S(  12) 

T  (  A  )  * 

S<  13) 

C  SOLVE  FOR  PARAMETERS. 

CALL  SOLVE 

IE  I  IERR  ,£U.  I )  RETURN 
A  *  AH 
8  »  8H 
C  *  CH 
U  *  OH 
GO  TO  30 
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20  Hli, i)  -  SI5) 

H ( 1,2|  *  5(6) 

H(l, 3)  *  SIT) 

H(2,l)  *  H (  1,2) 

H( 2 ♦ 2)  *  5(8) 

H(2,3i  *  S ( 9 ) 

HI  3,  li  *  H(  1,3) 

H( 3,2)  -  M{ 2 , 3 ) 

HI  3  » 31  *  SUO) 

Til)  *  SU1)  -  (A  *  S(  2 )  ) 

T*2)  »  S(  12)  -  (A  *  SI  3}  ) 

1(3)  »  S  (  1 3 )  ~  (A  *  SJA  )) 

C  SOL,:  FOR  PARAMETERS,. 

CALL  SOLVE 

If  I IERR  .EC.  U  RETURN 
a  a  AH 
C  *  BH 
D  »  CH 

C  COMPUTED  VALUES  OF  Y  AND  Y  RESIDUALS 

30  CO  AO  I  *  l,  N 

YCUS  »  A  ♦  IB  *  XIII)]  *  <c  *  X2U)J  ♦  ID  *  X3II ) ) 
YOEVI  I)  *  Yin  -  YCI  n 
AO  CONTINUE 
RETURN 
END 


ooonooon 
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IIBFTC  PARA 

SUBROUTINE  PARA 
COMMON 

C 

c 

COMMON 
COMMON 

C 
C 

COMMON 

C 

dimension 

C 

EQUIVALENCE 

C 
C 

c 

SUBROUTINE  FOR  DETERMINING  LEAST-SQUARES  SOLUTIONS  OE  PARAMETERS 
FOR  PARABOLIC  EQUATIONS  CF  FORM 

Y  *  A  ♦  (B*X)  ♦  1 C*f  X**2 ) I • 


A,  Al,  8,  81.  C,  Cl.  D.  Dlt  AH.  BH,  CH,  OH.  AN.  ANI, 
AN2.  OF  I .  OF2»  OFT.  FVALUE »  CO.  CV»  R,  POEVM.  SEY, 
XV.  YV,  YOEVSQ.  EA 

IA.  101 S  M »  I  ERR »  i F  V .  JSOIVE,  ISCI*  N,  NIV,  NOTE.  NP 
H14.4),  T 1 4 ) .  SDEV14J,  VMEAN14),  S130),  VC201.6), 
PDEVI20G)  .  X 1L  1 2'10 )  »  X2LI200).  X3U200I,  YU200I. 
YGI200),  Y0EVI200) 

! SC  14) .  JX ( 6 ) .  f MT { 18) .  TITL6116).  INO,  I  PAGE .  IEQ. 

I  CRD,  NIC*  NRUN 

1011201).  102(201),  Y1201).  XK201  ),  X21201), 
X3120I),  Fit) 

(Fin,  a$,  noun,  vii, u).  (102(1),  vu»2>), 

( Y  ( 1 )  ,  V  ( 1 , 3)  )  ,  1XUU,  V(  1 , 4 }  )  ,  1X2(1),  V  f  1 , 5 ))  » 

(  X3  ( 1 )  .  V(  1  *  6 )  )  ,  UX.l),  Jl).  (  JX  (  2  l  ,  J2)  , 

1JXI3),  J3),  I  JX  {  4  )  ,  J4 )  »  (  J  X  1  5  )  ,  J5),  UX16),  J6) 


C 


C 


WHERE 

A  MAY  BE  SPECIFIED 

IF  (  IA 

•EQ.  n  GO 

TO  20 

hi  i.n 

*  AN 

HU,  21 

-  SI 2) 

H( 1,3) 

-  SI  5) 

M2,n 

*  HU. 2) 

H12.2I 

«  S(  5) 

HI  2, 3) 

«  SUB) 

HI  3, 1) 

*  Hll. 3) 

H(  3,2) 

*  HI  2,3) 

H(  3,3) 

*  S(H) 

T  (  I )  » 

SID 

T(  2)  » 

sun 

TUI  m 

S 1 20 1 

SOLVE  l 

FOR  PARAMETERS, 

CALL  SOLVE 

IF  (  ICRR  «EQ.  1) 

RETURN 

A  *  AH 

B  *  BH 

C  «  CM 

GO  TO 

30 

20 

Hii.n 

*  SI5I 

HI  1,21 

«  SUB) 

H12.ll 

•  HU, 2) 

HI  2, 2) 

v  SI  19) 

(III  - 

sun  -  (a 

*  S( 2 ) ) 

¥121  « 

SI  201  -  IA 

*  sun 

SOLVE 

FOR  PARAMETERS. 

CALL  SOLVE 

IF  HERR  ,£Q*  U 

RETURN 

fl  *  AH 

C  *  8H 

:  r-; 


COKPUHi)  VAlULS  GF  Y  AND  Y  RESIDUALS 
K  UD  40  I  x  ! ,  m 

*C*J!  *  A  *  **<*'>  ♦  cc  *  (xiu; 

YDt  V  (  X  >  -  VH)  -  YC  {  I  ) 

40  tuNTINUf 

LAlCULAlf  CGORO 1  NATES  OF  VERTEX  POINT 
XV  =  -B/ I  2.0  *  t ) 

YV  -  A  ♦  ( B  *  XV)  +  (C  *  XV  *  *VJ 
Kf  T URN 


ENG 


r>  o  r>  o  o  o  o 
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SlflfTf,  POWR 

iiJBROUTf  N£ 
COMMON 

c 

COMMON 

COMMON 

f 

W 

t 

COMMON 

C 

OS  MENS  ION 

C 

EQUIVALENCE 


C 

c 


POWR 

A  ,  Alt  B,  Bl,  C,  Cl,  0,  01,  AH,  BH,  CH,  DH,  AH,  AMI, 
A N2,  OF  1 ,  OF 2,  OFT,  F VALUE ,  CO,  CV,  ft,  PDEVM,  S£V, 
XV,  YV,  YDEVSQ,  E A 

IA,  iDfSK,  ! ERR ,  IFV,  ISOLVE.  1  SCI •  N,  N1V,  NOTE,  NP 
H(A»A),  T  (  A )  ,  SDE  V  (  A )  ,  VME AN|  A  )  ,  S  1  30  >  ,  VI20l,6), 
POEV(  200) ,  X  1 L  C  200  )  *  X2LI200),  X3L1200I,  YLI200I, 

YC ( 200) ,  YQE V ( 200 ) 

I  SC ( A ) ,  JXife),  FMT ( 18  J ,  TITLEI16),  1MD,  IPAGE,  IEQ, 

1 ORG,  NID,  NRUN 

101 (201 ) ,  102(201),  Y ( 20 1 ) ,  X l ( 201 ) ,  X2(201), 

X3( 20 1 ) ,  F  C 1 ) 

(fin.  A),  iioim,  vti,m,  1102m,  m,2n, 

(Y(l),  V ( 1 , 3  ) )  ,  1X1(11,  V( 1  ,A ) I  ,  |X2**5,  V ( 1 , 5 ) ) , 
IX3U5*  ni,  6)),  (JX(l),  Jill  (  J  X I  2  )  *  J2  )  , 

(JX(3>,  J3  )  ,  (  J  X  l  A )  i  A )  ,  (JX(5),  J5) ,  (JXife),  Jfe) 


SUBROUTINE  FOR  OETERMINING  LEAST-SQUARES  SOLUTIONS  OF  PARAMETERS 
FOR  POWER  EQUATIONS  OF  FORM 

Y  *  A  ♦  ( X 1**0 )  *  (X2**C)  *  ( X 3**0 ) 


C 

c 


SET  SUBROUTINE  INDICATOR  TO  2. 

10  1N0  -  2 

DO  20  I  *  1,  N 

SET  ALL  VALUES  OF  EITHER  THIRD  <X3)  OR  SECOND  AND  THIRO  (X2,  X3) 
INDEPENDENT  VARIABLES  TO  l  IF  NOT  BEING  CONSIDERED  FOR  THIS  RUN. 
IF  (NIV  .LT.  3)  X 3 ( I  I  *  1.0 
IF  (NIV  .LT.  2)  X 2 4  I  )  *=  1.0 
X2LU)  *  ALOG  (X21I)) 

X3L 1(1  *  ALOG ( X  3 ( I  )  ) 


H( 1,2) 

HU, 2) 

♦ 

XILII) 

H( 1,3) 

« 

HU, 3) 

* 

X2LI1  ) 

HI l,fe) 

ac 

H(l, A) 

*■ 

X3LII  ) 

Hi  2,25 

Hi 

HI  2 , 2 ) 

♦ 

IX1LI  I  ) 

♦ 

XIU  I 

H(2, 3) 

4Y- 

H  (  2 , 3  ) 

+ 

(XILII) 

X2L  1  I 

HI, 2. A) 

Vi 

HI2, A) 

♦ 

(XIU  I  ) 

• 

X3LU 

H  (  3  *  3  T 

iz 

HI  3, 3) 

♦ 

(X2U  I) 

• 

X2LU 

H(  3  « A  ) 

M 

HI  3, A) 

♦ 

IX2LU  ) 

• 

X3LU 

HI A,AI 

ar 

MIA, A) 

*■ 

IX3U  II 

* 

X3U  I 

if  ( 2  ?  *  M2)  ♦  (XIL(i)  ♦  YL(I)  > 

T  (  3  5  *  T  (  3 1  ♦  ( X2L (  1  )  ♦  YL  (  I)  ) 

TlAi  *  HA)  ♦  (  X  3L  (  I  )  *  YL  4  I  I  I 

20  CONTINUE 
HU,  1)  *  M 
Hi  2 , 1)  *  H5 1,21 
HI  J,U  ■  Hi  1,3) 

Ml  3,2)  -  Ml  2, 3) 

HI  A,  n  *  Hlh  A’/ 

H  (  A  ,  2 )  *  H  (  2 »  A I 
Hi  A, 3)  *  M( 3„A) 

TUI  -  S  <1  A  * 

FIRST,  DETERMINE  LEAST-SQUARES  SOLUTIONS  OF  PARAMTERS  Al,  Bl,  Cl, 
01  FOR  LOGARITHMIC  FORM 


o  o 
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LN<  Y )  *  LN(Al)  ♦  ( B I  *  )  ♦  (Cl  *  I.  N  i  X  2 )  >  *  (01  *  LNIX3H 

CALL  SOLVE 

IF  (  IERR  .EC.  I >  RETURN 
IF  (AH  .GT.  8S. )  GO  TO  30 
A1  -  E Y  P ( AH ) 

Bl  «  BH 
Cl  *  CH 
01  -  OH 

C  DETERMINE  lEAST-SQUARES  SOLUUONS  OF  PARAMETERS  A*  5,  C»  0 

C  FOR  EXPONENTIAL  FORM. 

CALL  ITER 
RETURN 

r  '-RITE  ERROR  ME S5AGE. 

30  WRITE  ( 6 »  *G»  *h 

AO  FORMAT  (IHO//,  IOX,  86HTHE  EXPONENT  AH  IS  GREATER  THAN  88.  i H 

L  SUBROUTINE  PGWR.  THIS  RUN  HAS  BEEN  TERMINATED  '//  IOX, 

C  SHAH  *  ,  F15.7  > 

IERR  *  I 

RETURN 

END 


i 

# 

=s 

3 


,  i 

i 

i 


o  o  r>  o  o  r>  o 
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6IBFTC  A$VH 

SUA&OMTINE  A$¥M 
COW WON 

c 
c 

COMMON 
COMMON 

C 
C 

COMMON 

C 

DIMENSION 

C 

EQUIVALENCE 

C 
C 
C 

DIMENSION 


At  AI»  »,  81,  C,  Cl*  U »  Oi ,  AH,  RH,  CH,  DH,  AN,  ANi, 
AN2 ,  DF1,  DF2,  Of*,  FVAIOE,  CD,  CV*  R.  PDEVM*  SEY* 
XV,  TV,  YSEVSQ,  EA 

I A ,  IPTSX*  rfSR,  I F  V ,  ISOLDE*  ISCI,  N*  *»IV,  NOTE,  HP 
H(4,4}»  T  (  4 )  ,  SDEVK),  VH6AN14I,  $»39),  V(201*6), 
PDEVI20QI*  X  IL  (  200  f ,  X2l(200i,  X3LI200I,  v|,?00), 
YC120Q),  VD6VI200I 

I  SC  ?  4)  •  JXS65,  FMT(lfiJ*  TITUdA),  IND,  1  PAGE  •  !EQ* 
iORO*  NID,  HR UN 

101(2011*  1021201),  V(20D,  Ill  201 1  *  X2<201), 

X3(20i I,  F(  U 

(  F  (  1  i  r  A).  IlfiHl).  (102(11*  VI  1  »2 1 S  , 

tvm,  vii,s5J,  ixim,  vu*«n*  (X2U>,  vu*iD* 

mm,  v(i*6»i,  {j im*  ill,  uxi2i»  J2i» 
tJXO),  j> )  t  (JX(4),  J4 ) »  (JX(S>,  J5)  *  (JX(U,  Jfc) 
SUM  16),  bSAVEOI 


SUBROUTINE  FOR  DETERMINING  LEAS T- SQUARE S  SOLUTIONS  OF  PARAMETERS 
FOR  ASYMPTOT I C- POWER  EQUATIONS  OF  FORM 


Y  *  (A  ♦  t  X**B » J  ♦  C 


SET  INITIAL  VALUE  OF  0  TO  -4.01. 

10  Bi  »  -4.01 
C  SIT  FIRST- ITERATION  DESIGNATOR  TO  1. 

ITERS  -  1 

20  DO  190  K  «  l*  160 
A3T0RE  -  0* 

RSTOWE  «  0. 

CSTORE  •  0. 

C  SET  B  INCREMENT  INITIALLY  TO  0.05 

JO  DM  «  O.OS 

C  STEP  INITIAL  0  VALUE  BY  DB  INCREMENT. 

Bl  -  Bl  ♦  DB 
IF  (K  . NE -  01)  GO  TO  40 
ITERS  •  l 
Bl  -  B.01 
40  8  •  Bl 
$0  DS  60  I  *  I*  6 
SUM  1  II  •  •• 

60  CONTINUE 
70  DB  80  I  -  l*  N 
XP  •  XKII66B 
XPSQ  •  XP  •  XP 
XP  1  ■  XP  4  XIIU  ) 

SUM(  1)  *  SUW( 1)  ♦  XP 

SUM ( 2 1  -  JUM( 2 1  ♦  XPSQ 

SUM13I  •  SUMO)  ♦  XPl 

SUM  (  4)  -  SUM  (  4  )  *  (XPSQ  XILID) 

SUM(S)  -  SUM(  5 )  *  (Y(  I!  •  XPl) 

SUM(6l  •  SUM ( 6 )  ♦  ( ( v | | |  -  VMEAN(l))  •  XP) 
80  CONTINUE 

A  •  SUM(6)/( SUM! 25  -  (lUMfll  *  SUM) 1 ) / AN ) ) 
IF  (1A  .EQ.  U>  C  «  (S(l)  -  (A  •  SUM ( 1 ) ) ) / AN 
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G  *  SUM  (  5  i  -  (A  *  SUM!*,))  -  (C  *  SUM(3)> 

IF  ( K  .FQ.  1)  QSAVE(l)  =  SI21)  -  (A  *  SUM(6>> 

IF  «K  .LU.  160)  QSAVEI3)  =  S<21>  -  (A  *  SUM<6>> 

IF  (ITERS  .60.  2)  GO  TO  120 
IF  ( Gl  90,  170,  100 
90  M  -  -I 
GO  TO  110 
100  M  *  1 
110  ITERS  -  2 
GO  TO  190 

110  IF  ( M  .GT.  01  GO  TO  130 
IF  (Gl  150,  170,  UO 
HO  IF  (Gl  140,  170,  150 

no  b  *  B  -  (OB  *  0.5) 

GO  TO  160 

150  IF  (DB  .GT.  O.OA)  GO  TO  190 
B  *  B  «■  (OB  *  0. 5) 

160  OCA  «  A«S(A  -  •  < TORE ) 

OOB  *  ABS ( B  -  b STORE ) 

ODC  *  ABSIC  -  C  STORE  ) 

IF  (ODA  .LE.  1.0E-C8  .AND.  DOB  .LE.  1.0E-08  .AND.  ODC  .LE.1.0E-08) 
C  GO  TO  170 
ASTORE  =  A 
BS T ORE  *  B 
CS'ORE  «  C 
UB  *  OB  *  0.5 
GO  TO  50 

C  USE  NEW  VARIABLES  FOR  TEMPORARY  SOLUTION. 

170  AA  «  A 

BB  *  B 
U  «  C 

C  SUM  OF  SUUARES  OF  Y  RESIDUALS 

YDEVSO  *  S ( 2  1 )  -  (AA  •  SUM( 6  I  ) 

IF  ( I  SOL VE  .EU.  0)  GO  TO  180 

IF  (YQFVSy  .LE.  QSAVE ( 2 ) )  GO  TO  180 

AA  -  ASAVE 

BB  »  BSAVf 

CC  *  CSAVE 

YDEVSO  *  OS A VE ( 2 ) 

C  STORE  PARAMETER  -/ALUES  AND  SUM  OF  SUUARES  OF  Y  RESIDUALS. 

180  ASAVE  *  AA 
BSAVE  -  BB 
CSAVE  »  CC 
USAVE (2)  •  YDEVSO 

C  St T  FIRST-! ration  DESIGNATOR  Tn  1. 

ITERS  -  l 

c  SET  soldi  ION  DESIGNATOR  TO  I. 

1  SOLVE  *  I 
190  CONTINUE 

C  IF  A  'JNIOUE  SOLUTION  FOR  p  DOES  NOT  FUST  (N  THE  SPECIFIED  RANGf, 

(.  R  R  I N 1  A  MESSAGE  RELATING  TO  THAT  FACT. 

IF  tlfOlVE  .FC.  1  .AND.  YOEVSO  .LT.  uSAvEIU  .AND. 

C  YDEVSO  .11.  OS^VEim  GO  TO  210 
kRITE  16,  200) 

VOO  FORMAT  (  1  HO/ / ,  9X,  75HNO  SOLUTION  HAS  BEEN  FOUND  FDR  TH'S  PROBLEM 
C  IN  THE  RANGE  0»  -A  TO  ♦<.  FOR  B.  ) 
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C 


c 


^10 


22  0 


2  ERR  •  1 
RETURN 

RESTORE  VARIABLES  TO  ORIGINAL 

A  «  AA 
B  «  B8 

C  -  CC 

DOL«c*”“  °'  ’  ‘ND  r  "ESIOU.tS 

YCt  !>  *  IR  •  DO(I)MB)  I  r 
TOE V(  I)  «  V(||  -  YC ( I  ) 

CONTINUE 

RETURN 

ENO 
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*  I  B  F  T  C  EXPO 

SUBROUTINE  expo 

COMMON  A,  Al,  R,  HI,  C,  Ci,  0,  01,  AH,  BH ,  CH,  DH,  AN,  AN  l  , 

c  AN? »  DM,  OF?,  T,  TVALUE,  CD,  CV,  R ,  PDEVM,  SEY, 

c  XV,  YV,  YDtVSC,  EA 

COMMON  IA,  (DISK,  !  f  R  R ,  JFv,  ISCLVE,  ISCL,  N,  NIV,  NOTE,  NP 

COMMON  H(4,4),  T(4),  S0EV(4),  VME  AN( 4 )  ,  S(30),  V(201,6I, 

C  PDE  V  <  200) ,  X  l  L  (2001  ,  X?L(?OOI,  X3K200),  YL(200), 

C  YC ( 2C0 ) ,  YDEVI200J 

COMMON  l  $C  (  4 )  ,  JX  (  6  )  ,  FMTIlh),  T  1  T  1.  E  (  1  6  )  ,  1N0,  !  P  AGE  ,  1  E  Q  , 

C  IORU,  NIC,  NRUN 

DIMENSION  101(201  1  ,  102(2013,  Y(20l),  XU201),  X2I201), 

C  X  3 ( 2  C i ) »  Ell) 

EQUIVALENCE  (F(U,  A),  ((01(13,  w  (  1,13),  (102(1),  V(l ,?)), 

C  (Y(l),  V(),))l,  ( X  1  (  i  >  ,  V (  1,4)  ) ,  ( X2 (  1 ) ,  V(l,5)l, 

(  (  X  3  (  1  >  ,  VI  1,6)  )  ,  <  J  x  (  1  )  ,  J  1)  ,  (  J  X  (  ?  )  ,  J  2  )  , 

C  ( J  X ( 3 )  ,  J  3 )  ,  (  J  X ( 4 )  ,  J  4  t ,  (JX(S),  J5),  (  J  X  (  b  )  ,  J  6 ) 

C 

C  SUBROUTINE  f  OR  DETERMINING  L E A S 1 -SQUARE S  SOLUTIONS  OE  PARAMETERS 

C  FOR  EXPONENTIAL  EQUATIONS  CE  FORM 

C 

C  Y  •  (  X P ( A  ♦  (B  *  Xj| 

C 

C  (1RST,  OBTAIN  L E AS  T-SQUARE S  SOLUTIONS  Cc  PARAMETERS  Al,  Bl  FOR 

C  SI  M l-LCGAR I THMIC  FORM 

C 

C  L  N 1  Y  )  *  A 1  ♦  <  B 1  *  X  > 

C 

C  SET  SUBROUT  INF  INCICATOm  TO  ?. 

1C  INU  -  2 

C  CALCULATE  DENOMINATOR  OF  Al  TERM. 

UENOM  .  (AN  *  S(S))  -  (  S  (  2 )  •  S  (/••)) 

i.  CHECK  FOR  7(R0  Of  NOMINATORS. 

IE  10EN0M  .£(..  0.  .OR.  S ( 2  3  .EQ.  0.)  GO  TO  20 
(.  CALCULATE  NUM|RATO«  OE  Al  TfcKM. 

ANUM  ,  (SIS)  *  S'14J)  -  ( S ( 2 )  *  M!6)» 

C  CAL  cut  AM  Al  TERM. 

Al  *  ANUM/UENOM 
CAl CUE  ATE  Bl  TERM. 

Bl  -  (  S (  l*  '  -  (AN  •  A  1 )  )  /  S ( 2  ) 

C  DETERMINE  EE AST-SCUARES  SCU  T  IONS  Or  PARAMFTERS  A,  R  FOR 

EXPONENTIAL  FORM 
(  ALL  I  T  F  0 
RET  URN 

C  -RITE  FRROK  M(SSAG£. 

2C  »R  I  T  f  (6,  30  1  OF  NON,  SI  2),  AN 

3 C  FORMAT  (  1H0// , 10* , 11  1HA  SOLUTION  CANNUT  BE  OBTAINED  BECAUSE  A  ZERO 
(  DENOMINATOR  EXISTS  IN  CNF  CR  MORE  OF  THE  TERMS  (SUBROUTINE  EXPO). 
C  III  10X,  BHOENOM  »  ,  FIS. 7  /  10x,  «HS(2)  «  ,  FIS. 7  /  10X, 

(  KHAN  »  ,  ►  1 S .  7  3 
(ERR  *  l 
ME  TURN 
END 
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*  I  8F  re  SOLVE 

SUBROUTINE  SOLVE 

COMMON  A,  A 1 


COMMON 

COMMON 


COMMON 
DIMENSION 
CIMENS ION 

equivalence 


A,  Al,  8,  81,  C,  Cl.  0,  01,  AH,  8H,  CM,  OH,  AN,  AN l , 
AN2,  DEI,  OF?,  OFT,  f VALUE ,  CO,  CV,  R,  PDEVH,  SEV. 
XV,  Y V,  YDEVSU,  EA 

IA,  IOISK,  IERR,  iEV,  ISOLVE,  I  SC  A ,  N,  NlV,  NOTE,  NP 
H ( A , A )  ,  T { A  )  ,  SOEV(A),  VMEAN(A),  S  (  30  )  ,  VI20U6), 

POE  V  1 200  ) ,  X  l  L  (  2  0(J )  ,  X2K  200M  X3LI200),  YL(200), 
YCI2C0),  YDEVI200) 

ISC(A),  J  X  1  6 1  ,  F  m  T  (  1 0  )  ,  1  1  f  L  t  (  16),  IND,  IPAGf,  IEQ, 
1080,  NIL,  NRUN 

101  (201),  102(201  ),  YI201),  XM201),  X2I201W 
X  3 ( 201  )  ,  F  (  1  ) 

U( 12)  ,  *120 

!  F  (  1  )  ,  A),  (I  OKI),  V(l,ll),  (102(1),  V  (  1 , 2  S  )  , 

(Y(l),  V ( 1 , 3 ) )  ,  (Xl(l),  V (  i  ,  A  )  )  ,  (X2(l»,  V(1,S)», 

(  X  1  (  1  (  ,  V  (  1 , 6  )  )  ,  (jxm,  Jl),  (  J  X  (  2  ) ,  J  2  )  , 

! J  X (  1 )  ,  Jl),  ( J  X ( A  )  ,  J  A )  ,  (JX(S),  JS>,  ( J  X ( 6 ) ,  J  6  ) 


SUBROUTINE  FOR  SOLVING  SIMULTANEOUS  EQUATIONS 

GO  TO  (20,  10,  AO,  SO),  NP 
2  0  If  <H<  l  ,  I  )  .EG.  0.  !  ,0  TO  ftC 

AH  =  T  C  1  ) / H (  1 ,  l  > 

*  E TURN 

1C  OENCM  *  (Mil,  l ) /HI  1,2)5  -  (H  (  2  ,  1 ) /H( 2,2  )  ) 

IE  l DENOM  »  E  Q »  0.)  GU  TO  60 


ANu* 

■ 

(  T  (  I  t  /H  (  1  , 

)  -  (  T  (  2  )  /  h  (  2 

,2)5 

AH  * 

ANUM/'DENOM 

8H  - 

(  T  (  l  I  -  (AH 

0 

H(l;l)j)/H(l  ■ 

2  » 

RE  TUS 

N 

AO 

t(  l  ) 

M  {  M  l  )  • 

*■ 

X 

t 

X 

1  ) 

♦ 

H  ( 

M3)) 

U!  2  t 

i  H  |  1  ,  2  )  * 

hi  2 ,  1 )  )  -  ( H |  2  , 

2) 

* 

HI  i  ,  n  ) 

Oil) 

(Til)  *  H  (  2  » 

1 )  )  -  (  T  (  2  )  ♦ 

H(  1 

,  1) 

) 

GI  *  1 

DM  2,  l  )  ♦ 

H  ( 

1,1)1  -  t  H (  1 , 

l  ) 

0 

M (  2 ,  1)  ) 

U<  S  > 

(  H!  (  2 , 2  )  ♦ 

H(  3,  3)  1  -  ( H (  1  , 

2) 

0 

H  (  2  ,  1  )  ) 

U(  6  ) 

a 

(  T  !  2  )  *  H  ( 

1, 

1 ) )  -  (Til)  * 

H  ( 

) 

i 

,  1) 

i 

DENOM 

- 

(  U  (  l  )  /  G  (  2  )  ) 

-  (  u  (  A  )  /  u  1  S  » 

) 

IF  (DENOM  .EG.  0. 

) 

GO  TO  SC 

ANUM 

■ 

I U l  1 ) /UI 2  )  ) 

-  (U(6)/U(S)  ) 

AH  * 

ANUM/OENOM 

8H  « 

(U (  n  -  (AH 

• 

U  f  1  )  )  )  /  u  (  2  ) 

C  f  * 

<  T 

(  l  )  -  1  AH 

* 

H t  1  ,  1)  )  -  (  8H 

♦ 

H  (  I  , 

2  )  )  )  2H 

RETURN 

SO 

L(l) 

■ 

!“l 1 , 1 !  • 

H  i  2  .  A  )  )  -  (  H  1  2  , 

1  ) 

HI 

M  A  )  ) 

1(2) 

* 

(  H  I  1  ,  2  )  * 

H  ! 

2  ,  A  )  )  -  (  h  (  2  , 

2  ) 

Hi 

l  ,A  '  ) 

U(  1  1 

• 

0(1,1)  • 

H  ! 

2  ,  A  )  !  -  *  n  i  2  , 

1  ) 

m  i 

1,01 

D(  A) 

* 

(Til)  * 

H  ( 

2  .  A  )  t  (M2) 

H  ( 

1  ,  A  )  ) 

U(  s ) 

m 

!H(  1,  1!  * 

H  ( 

A  ,  A  >  )  -  (  H  (  1  , 

A  ) 

H 

3  ,  A  )  I 

u  ( * ; 

(HI  2,  1)  * 

! 

a,  A  M  -  0 1(2, 

A  ) 

H( 

3, A  I 

u(  n 

* 

(  H  1  1  ,  1  )  • 

- 1 

•»  ,  A  )  )  -  (  H  |  1  , 

A  ) 

Ml 

1  ,  A  )  ) 

<  i  n ) 

It 

nin  * 

"  1 

A  ,  A  J  )  -  |  T  (  A  ) 

HI 

3,  A  )  > 

c<  <n 

• 

(  H  (  i  .  2  1  • 

H  l 

•  ,  M  1  -  ( Hi  1  , 

A  I 

H( 

2  *»  M 

«  (  1C  ) 

* 

IM2.2)  * 

* 

l  •«.<•)  )  -  (  H  1  2 

,  A  I 

*  H  (  2  ,  A  )  ) 

u  m 

m 

M(  1,2)  * 

> 

!  A  ,  A  )  )  -  (  H  (  ( 

,  A  ) 

*  H  (  2  ,  A  )  ) 

U  12) 

r 

(M2!  • 

> 

»*  a  ,  A  )  )  -  MM 

) 

•  H  (  2  ,  A  )  ) 
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MU)  *  U(  7)  /UI  3) 

M(2)  a  UI 7 ) /U ( 6  I 
Ml  3)  =  U(  m/U(  71 

u ( 4 1  *  :ut4)/uiion  *  ( u ( 7 ) /u i ? ) ) 

M « 5  )  *  kill  -  I U I  8 ) /UI 4  1  ) 

M<6)  *  HI  3)  -  «U(  iOI  /*J<  €»>  1 
i*  4  7  )  ■  MI2)  *  I  UI 8) /UC 1 0) ) 

M(8)  «  Ml 1)  -  (U(6)/U<2) ) 

M(9)  *  M(3>  -  (UI  12  )  /U<  8  )  } 

-I 10}  *  UI  l )  /  U I  2  I 

Mil!  *  hill  -  tu(5)/u;  m 

h(17l  *  (U<  lll/UI  10)  >  -  M<2) 

HID)  -  (U<  5  >/U<  6)  J  *  iu(7)/ulicn 
hlH)  *  Will  -  I  UI  6 )  /Ut  2  I  ) 
hlh)  *  3)  -  <U(9)/Ul  51) 

Ml  165  *  M(4)  *  W(5)  ♦  W  <  6 ) 

Ml  17)  =  HI 7)  ♦  W I  8 )  *  H(9) 

DENOM  =  (Ml  10)  *  M  >  '  l )  *  M(12}}  -  (MI13)  *  W(l4)  •  HI15)) 

IE  ( DENCH  .FC.  0.}  GO  TC  60 
ANUM  *  Ml  16)  -  WU7) 

AH  *  ANUM/DENOM 

DENOM 1  *  lui?)/unn  -  (0165/0(71) 

IF  I  DENOM  l  .  EQ.  0.)  GC  TO  60 

mi  is)  =  ium/uni)  -  in 1 8 ) /ui 7 j ) 

Ml  191  *  (UI  1  )/U( 35  -  i  I  5  I /UI 7 ) * 

ANUM 1  »  hllRI  -  I  AH  *  Ml  195  1 
BH  =  ANUMl/OENOMl 

CH  =  |U<4)  -  (AH  *  Ul  1)  )  *  (BH  *  U(2))»/Um 

UH  -  (Til)  -  (AH  *  HI  1*15)  ~  1 BH  *  h(l,?l)  -  (CH  *  H I  i , 3 1 1 5 /H I  1 , 4 ) 
Kf TURN 

0  ERROR  MESSAGI 

60  MR»TE  (6,  70) 

70  FORMAT  (1M0//10X,  97HA  ZERO  DENOM!'  UOR  EXISTS  IN  THE  0  ALCUL AT !  ONS 
C  OF  SUBROUTINE  SOLVE.  THIS  RUN  HAS  BEEN  TERMINATED. 5 
GO  TO  (  120,  80,  100) ,  I UD 
80  MKSTE  (6,  90) 

90  FC.MAT  I IHO/  10X,  felHNOTE.  SUBROUTINE  SOLVE  HAS  LAST  CALLED  FROM 
(  SUBROUTINE  ROHR.  ) 

GO  TO  120 
100  MRITE  16,  110) 

110  «-URMAT  (  IHU/  10X,  61HN0IE.  SUBROUTINE  SOLVE  HAS  LAST  CALLED  FROM 
2  SUBROUTINE  I  TER.  ) 

120  JERR  *  I 
RETURN 
END 
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* ! BF  T C  ITER 

SUBROUTINE 

COMMON 

C 

C 

CONDON 

COMMON 

C 

C 

COMMON 


c 

DIMENSION 


C 

DIMENSION 

EQUIVALENCE 

C 

c 

c 


I  TER 

A,  Alf  8,  81,  C,  Cl,  0,  Dl,  AH,  8H,  CH,  DH,  AN,  ANl, 
AN2,  DF1,  DF2 ,  OFT,  FVALUE,  CO,  CV,  R,  PDEVM,  SEY, 

XV,  YV,  YDEVSQ,  EA 

IA,  IDISK,  I  ERR ,  I F  V ,  ISOLVfc,  ISC1,  N»  NIV,  NOTE,  NP 
HiA,4),  HA),  SUEV<4),  VHEAN(A),  SI30),  VI201,6), 
P0EVI200),  X  l  L ( 200 )  X2L1200),  *31(200),  YL1200), 

YC ( 200) ,  YDE V ( 200 ) 

I  SC ( A ) ,  JX ( 6  }  ,  FMT (  18),  T I  TIE ( 1 6 ) »  I  NO ,  IPAGE,  IEQ, 

I  ORO ,  NIO,  NRUN 

101(201),  102(201),  YI201),  X 1 ( 20 1 ) ,  X2I201), 

X  3 ( 20 l ) *  F <  1) 

r-P(A),  ATEMP  (  1 1  ) ,  8TEMPlll),CTEMPl  m,0T£MM(ll),0(ll) 
IF  SI),  A),  (101(1),  VU.D),  (102(1),  Vt  1,2)1, 
i  V  C  13,  V  (  1  ,  3  )  )  ,  mil),  V(  1 ,  A  )  )  ,  (X211),  V(l,5)5, 

(  X  3  (  1 ) ,  V  (  1 ,  A  )  1  ,  (jxm,  Jl),  (  JX  (  2  )  ,  J2 )  » 
lviX(3),  J3  )  ,  <  J  X  J  A  )  ,  JA),  IJXC5),  JS),  IJX16),  J6) 


C 

C  SUBROUTINE  FOR  DETERMINING  LEAST-SQUARES  SOLUTIONS  OF  PARAMETERS 

C  F8R  NON-LINEAR  EQUATIONS  WHERE  AN  ITERATIVE  PROCEDURE  IS 

C  REQUIRED 

C 

C  STEP  SUBROUTINE  INDICATOR  iY  l. 

1C  ! NO  *  i NO  *  1 

C  SET  INITIAL  GUESSES  TO  LOGARITHMIC  SOLUTIONS 

A  *  A1 
B  M  01 
C  m  Cl 
D  »  Dl 

20  CO  200  L  *  1,  50 
C  CLEAR  H  AWO  T  MATRICES^ 

CB  30  !  *  39,  58 

FID  n  0. 

30  CONTINUE 

If  I IEQ  .IQ.  3)  GO  TO  AO 

C  CHECK  MAGNITUDE  OF  A  FOL  EXPONENTIAL  CASE* 

IF  (A  .GT.  88.)  GO  TO  230 
EA  *  EXP ( A ) 

AO  CO  100  I  *  1,  N 

IF  ( IEQ  .EG.  3)  GO  TO  50 
El  *  A  ♦  IB  ♦  XI  (I )  ) 

E2  *  2.0  A  El 
E J  *  X  1L  I  1  )  ♦  E2 
EA  *  12.0  *  XlL(ll)  >  E2 
C  CHiCK  MAGNITUDES  OF  EXPONENTS. 

IF  I E2  .GT.  88.  -.OR.  E3  .GT.  68.  .OR.  6A  .GT.  88.)  GO  TO  2)0 
C  CON PUT EO  V  VALUES 

YCII)  *  EXP1EI ) 

C  DIFFERENTIALS  OF  Y  FUNCTION  WITH  RESPECT  TO  PARAMETERS  A  ANO  8 

FP1 1)  «  YCII) 

FP ( 2 )  *  YCII)  *  X 1 ( | ) 


GO  TO  60 


C  CHECK  MAGNITUDE  OF  TWICE  THE  PRODUCT  OF  EACH  PARAMETER  IB,  C,  U> 

C  TIMES  LOGARITHM  OF  VALUE  OF  INDEPENDENT  VARIABLE  FOR  WHICH 

C  PARAMETER  IS  EXPONENT.'  IF  GREATER  THAN  88.,  PRINT  ERROR  MESSAGE. 
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50  FAS  *  2 oO  *  fi  ♦  X1LU ) 

FAC  =  2.0  *  C  •  X2L I  I  ) 

FAS  »  2-0  *  D  *  XHI1J 

IF  (FAB  «6T.  8®.  »QR.  FAC  .GT.  88.  .OR.  FAD  ,GT .  88.)  GO  TO  230 
C  DIFFERENTIAL  OF  Y  FUNCTION  WITH  RESPECT  TO  PARAMETER  A 

FPUJ  *  (Xl(!)**h)  *  ( X2(  5 )  **C  i  *  (X3U!**D» 

Ycm  «  a  ♦  ppm  i 

C  DIFFERENTIALS  OF  Y  FUNCTION  WITH  RESPECT  TO  PARAMETERS  8*  Ct  AMD  0  i 

FP( 2 )  *  YC(I)  A  XiLI!)  I 

FPU)  -  YCI  II  *  X2L(  1} 

FP  f  4 )  *  YCUI  *  X3LCI) 

C  Y  RESI0UAL5 

6C  YDfVU)  *  vm  -  YCU)  j 

IF  (NOTE  .EC.  1J  GO  TO  IOC 
IF  <  i SOLVE  -EG.  0)  GO  VO  70 
YD  *  (2.0  *  YCU  II  -  Y(  I  ) 

IF  (YO  .LE.  0.1  NOTE  *  I 

GO  TO  100 

C  CALCULATE  H  ANO  T  MATRICES. 

70  00  90  II  «  1,  ♦ 

DO  80  JJ  *  1,  4 

Hlll.JJ)  *  HI  1 1 « JJ)  ♦  IFPUI)  *  FPIJJM 
80  CONTINUE 

Tciii  *  run  ♦  iyoevui  •  fpuim 

9C  CONTINUE 
100  CONTINUE 

C  IF  A  SOLUTION  HAS  BEEN  OBTAINED  USOLVE  *  II  ♦  STOP  ITERATION  AND 

C  RETURN. 

IF  I  I  SOLVE  .EG.  1)  RETURN 

C  SOLVE  FOR  CORRECTIONS  TO  PREVIOUS  SOLUTIONS. 

CALL  SOL VI 

IF  HERR  .EQ.  1)  RETURN 
DA  *  AH 
08  *  BH 
DC  *  CH 
CD  “  OH 
C 

C  FIND  WHICH  FRACTIONAL  PART  OF  CORRECTION  TERMSt  WHEN  ADDED  TO 

C  PARAMETER  VALUESt  GIVES  LOWEST  SUM  OF  SQUARE*  OF  V  RESIDUALS. 

C  * 

TEMP  *  i.O 

HC  TEMP  *  0.1  ♦  TEMP 

00  120  J  *  ly  11  ] 

FI  «  TEMP  ♦  FLOAT ( J  -  1) 

ATEMP( J I  *  A  ♦  (DA  *  FI )  ] 

BTSMP( J I  *  8  ♦  (08  *  F| I  I 

C  TEMPI J I  •  C  ♦  (DC  •  Fi I  I 

0TEMP( j)  »  0  ♦  (DO  *  FI  I  | 

0 ( J 1  •  0.  J 

l’C  CONTINUE  i 

00  160  J  *  I*  11 
DO  150  I  *  I,  N 
IF  < I EQ  .EQ.  5)  GO  .0  130 

YflMP«  ATEMP(  J»M  X 1 4  l)**8TiMP(  J))*(X2(  II**CTENP(  J)  I* 

C  (X3(I)**0TEMP(JI) 

GO  TO  140 


no  o  noon 
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i  30  YTfMP  »  EXPI ATENPIJ)  ♦  (BTfMPU)  *  Xim>) 
U()  YD  I F  ■  Y  { I  I  -  YTEMP 

fllj)  *  S I J  ?  ♦  ( VOIF  *  YOIF) 

15C  CONTINUE 
IfcO  CONTINUE 

YOfVSC  »  QU) 

L«  *  1 

90  1  TO  J  *  2*  11 

If  IVDEVSG  .IE.  OIJI)  GO  TO  170 

c#  **  J 

YDS VSQ  *  6IJ) 

170  CONTINUE 


IF 

(LM  .GT.  U  GO  TO  ICO 

IF 

I  ABS(DA*10»0*TEHP)  JGT. 

1.0E-08) 

GO 

TO 

no 

IF 

I  ABS ! 08*10. O^TEMP )  JGT. 

i.OE-08) 

GO 

TO 

no 

IF 

ABS(DCFIO.O«TEMP)  JGT, 

l.QE-08) 

GO 

TO 

no 

IF 

( AiS(DB*iO.O*TE«P)  JGT. 

1.0E-08) 

GO 

TO 

no 

180  ODA 

»  A9SIA  -  ATEMPILMH 

DOB 

»  ASS  IB  -  BTEMPRMI) 

QOC 

»  ABSIC  -  CTEMP(LHI) 

DOB 

*  ABS 1 D  -  DTEMPILM)  ) 

C  UP** TE  THf  VALUE  OF  EACH  PARAMETER. 

A  n  ATEMPILM) 

8  *  BTEMPRW) 

C  n  CTEMPILM) 

0  «  OTENPILMI 

IF  A  SO'UTJON  IS  OBTAI*  .D«  SET  SOLUTION  DESIGNATOR,  ISOLVC,  TO  1. 
(A  SOLUTION  IS  ASSUMED  WHEN  THE  CHANGE  IN  THE  VALUE  OF  EACH 
P»«AMETER  FROM  ONE  ITERATION  TO  THE  NEXT  BECOMES  EQUAL  TO,  OR 
LESS  THAN,  |0** (  —  * )  1 J 


190  Iff 

(ODA 

.0T» 

1 • OE-ifl ) 

GO 

TQ 

200 

IF 

( ODB 

•  §T  . 

l.CE-BB) 

GO 

fO 

200 

IF 

( ODC 

.GT. 

l.OE-98) 

GO 

TB 

200 

IF 

(ODO 

.LE. 

liOE-08) 

1  SOLVE 

»  i 

200  continue 


ERROR  MESSAGES 
210  WRITE  16,  220)  Al,  81 

220  FIRM AT  UNO//  10X,  96HN0  SOLUTION  HAS  BEEN  OBTAINED  FOR  THIS 

C  RUN  AFTER  90  ITERATIONS.  THIS  RUN  HAS  BEEN  TERMINATED.  /// 

C  10* ,  21HL0GARI THN1C  SOLUTIONS  /  10X,  5H  A  *  ,  F16.*  /  |0X, 

C  SH  8  •  ,  FI6.5  ) 

IF  (NIV  .EQ.  21  WRITE  16,  222)  Cl 
222  FORMAT  < IN  ,  9X,  9H  C  «  ,  F16.5) 

IF  iN?v  .|Q.  3)  WRITE  16,  226)  Cl,  Di 
226  FORMAT  UM  ,  «,  5H  C  m  ,  F16.5  /  IOX,  5H  D  »  ,  F14.9) 

GB  TO  290 

230  WRITE  16,  260)  A,  8,  C*  D 

260  FORMAT  (  IMO//  IOX,  66H0VcSFL0MS  EXIST  IN  SUBROUTINE  ITER. 

C  THIS  RUN  HAS  BEEN  TERMINATED.  //✓  10W,  6HA  -  ,  FI6.5  /  IOX, 

C  6HB  ■  »  f  16.9) 

IF  (NIV  .EQ.  2)  WRITE  (6,  262)  C 
262  FORMAT  <  IN  ,  9X,  6HC  •  ,  FI6.5) 

IF  (NIV  .iQ.  ))  WRITE  16,  266)  C,  D 
266  FORMAT  ( IH  ,  9X ,  6HC  *  ,  F16.9  /  IOX,  6H0  « 


F 16,9) 
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?5C  I f  R R  =s  1 
Rtf  URN 
LNn 
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SIBFTC  STAi 

SUBROUTINE  ST AT 

COMMON  A,  A 1 9  8 ,  81,  C,  Cl,  0,  01,  AH,  8H,  CH,  OH,  AN,  A  X  l  » 

C  AN? ,  OF l *  OF 2 ,  DFT,  FVALJE,  CO,  CV,  R,  PDEVM,  SET, 

C  XV,  W,  YOEVSQ,  EA 

COMMON  I A ,  SOI  SR,  IERR,  IFV,  ISOLVE,  I  SCI »  *4,  N!V,  NOTE,  NP 

COMMON  H  f  4,4 ) *  T  C  4 ) ,  SDE V ( 4  I ,  VME AN( 4 ) ,  S(JO),  V(20l,6), 

C  POE  V  4  200 )  ,  X  it.  (  200) ,  X2U200),  X3L12Q0),  YL(2Q0>, 

C  YC i ZOO  I ,  Y0EV(200) 

COMMON  ISC(4T,  JX ( 6) ,  FMT  t 18 )  ,  TITLE! 16),  INO,  IPAGE,  IE3, 

C  lOROv  NIO,  NRUN 

DIMENSION  101(201),  102(201),  Y(201l,  XK  201  ),  X2(201), 

C  «S(20U,  F  (  1 ) 

EQUIVALENCE  (Fill,  A),  ((Old),  V(l,D),  (102(1),  Vd,2li, 

C  (Yd),  V{1,3)),  mm,  V  (  1 , 4 )  )  ,  (X2d>,  V  (  l ,  5  )  )  , 

C  5X3(1),  Vll.61),  (JXd),  Jl),  (  JX  (  2  ) ,  J2 )  , 

C  (  JX  (3  (  ,  J  3  )  ,  (JXd)  ,  J4 ) ,  (JX(5),  J5  )  ,  (  JX  (  6) ,  J6) 

C 

C  SUBROUTINE  FOR  CALCULATING  STATISTICS 

C 

10  VDEVSQ  «  0. 

DO  20  I  *  1,  N 
C  PERCENT  Y  DEVIATIONS 

POE  V  ( I  )  *  (100.  *  YDEVUd/Vd) 

C  SUM  OF  ABSOLUTE  PERCENT  Y  DEVIATIONS 

POE VM  -  POE VM  ♦  ABSIPOEVdd 
C  SUM  OF  SQUARES  OF  Y  RESIDUALS 

YDEVSQ  «  YOEVSQ  ♦  (YOEV(I)  ♦  VOEVdl) 

C  SUM  OF  COMPUTED  Y  VALUES 

SI25)  «  SI25)  ♦  rCiil 
20  CONTINUE 

C  MEAN  OF  COMPUTED  Y  VALJES 

VCMEAN  *  $(  25 ) /AN 

C  SUM  OF  SQUARES  OF  COMPUTED  V  VALUES  ABOUT  THEIR  MEAN 

DO  30  I  «  1,  N 

S(?fe)  •  S(26)  ♦  (!YCd)  -  YCMEAN)  **2 ) 

30  CONTINUE 

C  TOTAL  DEGREES  OF  FREEDOM 

DFT  -  AN1 

C  DEGREES  OF  FREEDOM  ABOUT  REGRESSION  CURVE 

OF l  •  N  -  NP 

C  DEGREES  OF  FREEDOM  DUF  TO  REGRESSION 

OF 2  «  OFT  -  OF  l 
C  F  VALUE 

IFV  •  0 

OENOM  «  YOEVSQ/DFl 

IF  (OENOM  .NE.  0.  .AND.  DF 2  .Nfc.  0.)  GO  TO  4C 
IFV  •  1 
GO  TO  50 

40  ANUM  •  S(2fe)/DF2 

FVALUE  *  ANUM/OENOM 
(F  (FVALUE  .GF.  t.0FK)8)  IFV  *  l 
C  MEAN  OF  ABSOLUTE  PERCENT  Y  DEVIATIONS 

SO  PDEVM  ■  PDEVM/AN 

C  STANDARD  ERROR  OF  ME  ESTIMATE:  OF  r  (  ^OJUSTtt:  F  04  SAMPLE  Sdt) 

SI  ¥  *  SQR  T  (  YOF  VSj  /  )F  l  ) 
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C  COEFFICIENT  OF  VARIATIO H  PERCENT 

CV  *  <100-  *  SEYI/VMEANC  l) 

C  COEFFICIENT  OF  DETERMINATION  lUNADJUSTED  FOR  SAMPLE  SIZE) 

CD  «  1.0  -  IVOEVSQ/SC2I)) 

C  COEFFICIENT  0*  CORRELATION  IUNAOJUSTEQ)  FOR  SAMPLE  SIZE) 

R  »  SQRTICO) 

RETURN 

END 
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UBFTC  OUT  I 

SUBROUT! HE  OUT l 


COMMON 

A,  Bf  81,  C  j  Cl,  0,  D) *  AH,  BH,  CH,  OH,  AN,  ANl » 

AN2,  DF1,  OF 2,  DFT,  (VALUE  *  CD,  CV,  A,  PDEVN,  SEY, 
XV,  YV,  Y9EVSQ,  EA 

COMMON 

IA,  IDISX,  If  RR*  IF  V,  I  SOLVE ,  I  SCI  *  9.  NIV,  NOTE,  NP 

COMMON 

H ( 4, 4 ) ,  TI4),  SDEVI4),  VWEANI4),  SIBO),  V(201,6), 
F0EVI200),  X1LI200),  X2i(200),  X3L(200),  YU200), 
YCI2D0),  YDCVI20G) 

COMMON 

I  SC  1 4) ,  JXU),  FMf(lB),  TITLE!  16),  INO,  IPAGS,  IEQ, 

I 0*0,  NIO,  NR UN 

DIMENSION 

101(201),  102(201),  Y( 201 ) ,  Xl(20t).  X2I291), 

XBi  201 ) ,  FUi 

EQUIVALENCE 

(F(l),  A),  (101(1),  V(l,lli,  (102(1),  V( 1 ,2 ) ) , 

(Yd),  VII, 3)1,  (Xlil),  VII, 41),  (X2U),  V(l,5)), 
1X3(1),  V(l,A)l,  ! JX ( 1 1 ,  Jl),  1 JXI 2),  J  2 )  , 

( JX ( 3) ,  Jl),  ( JX ( 4 1 ,  J4) ,  ( JX ( 5  ) ,  Jt),  ( JXI A) ,  JA) 

C  SUBROUTINE  FOR  PRINTING  SUMMARY  TABLE 
C 


10  WRITE  (A,  201 

20  FORMAT  ( IHO//  59X,  1 3HSUMNARV  TABLE  ///I 
HPl  •  *  I A 

IF  (IA  »EQ.  II  aRITE  ik,  25) 

25  FORMAT  ( 1H  .  12X,  1Q4HMQTE  —  THE  STATISTICS  CALCULATED  FOR  THIS 
C  RUM  ARE  MOT  COMF ARABLE  WITH  THOSE  FOR  UNSFECIFIED  V-lNlERCEPTS. 

C  f  f  i  ) 


30  WRITE  I  A*  401  Av  B 
40  FORMAT  ( IW  t 
C  BOX •  41  HA 
C  30X,  41HB 

GO  TO  190,  *0,  50,  TO),  NF1 
50  WRITE  (A,  AO?  C 
60  FORMAT  I iH  , 

C  29X ,  41HC 

IF  (IEO  «EQ.  ■>)  WRITE  IA,  A5)  XV,  YV 
65  FORMAT  I IH  / 

C  BOX,  41HX  COORDINATE  OF  VERTEX  FOIMT 

C  BOX,  4IHY  COORDINATE  OF  VERTEX  POINT 

GO  TO  90 

TO  WRITE  16,  80)  C,  U 
80  FORMAT  ( IH  , 

C  29X,  41HC 

C  BOX,  41 HO 

90  IF  I IEQ  .NE.  B  .AMO.  IEO  .ME.  5)  GO  TO  190 
WRITE  (6,  1001 

100  FORMAT  UH  /  BOX,  21HL0GARI  THMIC  SOLUTIONS 
WRITE  IA,  110)  Al,  B1 
110  FORMAT  (IH  , 

C  29X,  41H  A 

C  BOX,  41H  B 

GO  TO  (160,  1 AO,  120,  140),  NP1 
120  WRITE  (A,  130)  Cl 
1 30  FORMAT  UH  , 

C  29X ,  4 1 M  C 
GO  TO  160 


,  10X,  Mi.)  / 
,  IOX,  F 14. *  » 


,  IOX,  F14.9  ) 


,  191,  F14.5  / 
,  IOX,  F14.5  ) 


.  IOX,  F14.5  / 
,  IOX,  F 14. 5  I 


•  )0x,  FI 4. 5  / 

*  IOX,  F14.5  ) 


.  IOX,  F 14. 5  ) 


-1*0- 


10X, 

ID*, 


10X, 


140  write  (6,  iso)  ci,  oi 

ISO  FORMAT  ( IH  , 

C  2 9x,  4  l  H  C 
C  30 X,  41 H  0 

160  If  (IEQ  .EU.  5)  WRITE  (6-  1/0)  EA 
l 70  FORMAT  ( IH  / 

C  30 X ,  4 1 H  Y  INTERCEPT 
190  WRITE  (6,  200)  R,  CO,  SEY,  CV,  YOEVSQ,  '’DEVR 
200  FORMAT  ( IH  / 

C  30X,  4 l HCOEFF ! C I E  NT  OF  CORRELATION  IUNADJUsTEO)  ,  10X, 

:  30X ,  4HC0EFFIC1ENT  OF  DETERMINATION  IUNADJU  TED),  10X, 

C  30X ,  46HS  T  ANDAHO  ERROR  OF  THE  ESTIMATE  OF  <  (  f.  uj  JS  TED  )  * 

C  F14.S  / 

C  30X,  41HC0EFF ic IENT  OF  VARIATION  (PERCENT)  ,  IDX, 

C  30 X,  4 1 HSUM  OF  SQUARES  OF  Y  RESIDUALS  ,  10X, 

C  30X,  '♦IHMEAN  of  ABSOLUTE  PERCENT  Y  DEVIATIONS  ,  IDX, 
SF  (  I F  V  .  EQ.  0)  WRITt  (6,  210)  FVALUE 
210  FORMAT  (IH  / 

C  30  X ,  4 1 HF  VALUE  ,  IDX, 

IF  (  I F V  .EO.  I!  WnlTt  <6,  220! 

220  FORMAT  1 IH  / 

C  30 X,  41HF  VALUE  ,  16X, 

L  30HEQUAL  TO  OR  GREATER  THAN  10**P  ) 

WRITE  16,  .30)  OF l ,  0F2,  OFT 
230  FORMAT  ( IH  , 

l  2  9  X ,  41-  SCREES  OF  FREEDOM  ABOUT  REGRESSION  CJRVE,  IDX, 

C  30X,  4 1 Hu  t  OR  f  F  S  OF  FREEDOM  DUE  10  REGRESSION  ,  IDX, 

l  30  X ,  41HT0ML  DEGREES  OF  FREEDOM  ,  IDX, 

WRITE  (6,  2  3SI 

23S  FORMAT  (IH  JOX,  19HMEANS  OF  INPUT  DATA  ) 

‘'RITE  (  6,  2401  VMEAN(l),  VMt  AN!  2 ) 

240  F OR M  *  T  (  IH  , 

C  29X,  41H  Y  ,  IDX, 

C  30X,  4  If-:  XI  ,  IDX, 

oO  I  0  1  290,  2 SO,  270)  ,  Nl V 
2  SO  WRITE  (6,  260)  VMf  AN (  3 ) 

260  FORMAT  ( IH  , 

C  29X,  4 1 H  X  2  ,  ID*, 

00  ID  290 

2/0  WRITt  (6,  2  HO  )  v  ME  AN  (  3  )  ,  VMEAN(4I 
2  BO  FORMAT  (  IH  , 

C29X,  41MX2  ,  l OX, 

L  30X,  4 1 H  X  3  ,  IDX, 

290  WRIT!  16,  29S) 

2  9  S  FORMAT  (H  /  3  0  X ,  33HSTANOARU  DEVIATIONS  JF  INPUT  DATA) 
WRITE  (6,  300)  SUE  V ( 1 ) ,  S0EV(2) 

300  FORMAT  (IH  , 

C  29X ,  4 l H  Y  ,  IDX, 

C  3  0X ,  4  1H  XI  ,  I  OX, 

00  10  (  3SU «  310,  3)0) ,  Nl  V 

310  WRITE  (6,  320)  SDfcVl  3) 

120  F  JRMAT  (  IH  , 

!  2  9  X  ,  4  1  H  X  2  ,  1  D  X  , 

,0  TO  )  SO 


3)0  wRIU  (6,  340)  bOE  V (  3  > ,  SDE  V ( 4 ) 

3  4  0  f ORMAT  (  IH  , 


F14.S  / 
FJ4.S  ! 


F14.5  ) 


F14.5  / 
F 14. S  / 
5  X  , 

F 14. S  / 
F14.5  / 
F14.S  ) 


F14.5  ) 


F14.5  / 
F14.S  / 
F14.S  ) 


F14.S  / 
F14.S  ) 


F14.S  ) 


F14.S  / 
F  l  4. S  ) 


F14.S  / 
F14.S  ) 


F14.S  ) 


r\  o 


II8FTC  GU T 2 

SUBROuTlN E 
COMMON 


COMMON 

COMMON 


cui  2 

A  ,  Aly  By 
AN2,  DU, 
XV,  YV,  Y 


,  Dl,  AH,  0H,  CH,  OH,  AN,  AN  1  , 
ALUE,  CO,  CV-  R,  F'DE  VM t  SEY, 


HI,  C  i  C  !  ,  C),  Dl,  AH,  0H , 

AN2 »  DU,  0F2  »  OFT,  F  VALUE  ,  CO,  CV-  R,  F'DEVM*  SEY, 
XV,  YV,  YDFVSa,  EA 

(A,  (DISK,  IERR,  i  F  V ,  ISOLVE,  I  SC  1  •  N •  N1V,  NOTE.  NP 
H(A,A),  T(A),  SDEV(A),  VMEAN(M,  S  f  30 ) ,  VI201,6>, 
PDEVI200),  X1LI200),  X2L(200),  X1LI200),  YL(200), 
YCI2C01,  Y0EV12O0) 

(SClA),  J  X  (  6 }  •  F  H  T  (  IP),  TITLE!  16),  INC),  IPAGE,  1EQ, 

I CRC ,  N1D,  NRUN 

10  1(201  ),  102(  201),  Y ( 2  0 1 )  ,  X  11201  ),  X2(201), 

X  31 7CI  )  ,  F (  1) 

( F  (  1 )  ,  A),  (101(1),  Vll.m,  (102(1),  V  (  1,2)), 

(Y(l),  v<I,3)),  (XI  U),  V  (  1 ,  A  )  I  ,  (X2(l),  V  (  l ,  5  H  » 

(  X  3  (  1 )  ,  V<1,6»),  IJMII,  Jl),  (  J  X  1  2  ) ,  J  2), 

( JX ( 3 )  ,  J  3 ) ,  (JX(A),  JA),  ( JX ( S ) «  J5I ,  (JX'6),  J6) 


COMMON 


C I  HE  NS ICN 


EUU1 VALENCE 


SUBROUTINE 
Y  RtSIOUAl 


FOR  P  0  I  N I  I  No  INPUT  DATA,  CAL  CULATED 
S,  AND  PERCENT  Y  DEVIATIONS 


VALUES, 


PRINT  TITLE  ON  NEW  PAGE 
1C  WRITE  (fc,  2  0  ’  (TITLEU),  I  -  1,  16),  IPAGE 
20  FORMAT  <  1 H  i  /  1 0  x ,  16AA,  A l X  ,  5HPAGE  ,  12  / 
STEP  PAGE  NUMBER  By  i. 

IPAGE  '  IPAG(  ♦  1 

St  T  LINE  COUNT  tn  2FR0. 

LINF  S  -  C 
WRITE  (6,  30' 

30  FORMAT  (  1  HO /  A  9  X ,  UHCCMPUTFC  y  VALUES  AN 
GO  TO  (AC,  luC,  1  AC )  ,  N 1 V 
A (  WRITE  U,  SC) 


VALUES  AND  RESIDUALS  ) 


SC  FORMAT  (luO/  mvx,  TMPtRCiNT 


l 

1  s  X  , 

THY 

CALC.,  1  3  X ,  2 ( fcH  Y 

UO  90 

l  « 

l  ,  N 

F  I  1  - 

I  M 

F  1  2  3 

F 1  OAT (  1  >  /  S.C 

(F  (LINES 

.LI.  AO)  GO  (0  SS 

wR  (  TF 

(6. 

2  0  1  (  T  l  I  L  F  (  X  )  ,  T. 

1  PACE 

«  I  PAGE  ♦  1 

WRITE 

(  6  , 

p 

wR  I  TE 

(  fc, 

St  1 

II  NT  S 

*  C 

ss 

WRITE 

<6* 

60 !  ( V (  I  *  X )  ,  K  - 

(  c 

F  UR  MAT 

(  1H  ,  l  3X  ,  2  A A ,  S | SX 

t  1  NF  S 

*  l 

j  N  i  s  ,  1 

If  (  f  11  . 

1  ..  M2)  WR  I  Tfc  !  f-  , 

SHL  ABLL  , 
I X )  /  ) 


1  7  X ,  1  HY  t  l  7  X  ,  2  H  X  i  , 


-  1,  L  6 )  ,  (PA  ,c 


1  ,  <•  )  ,  YC  (  I  )  ,  Y  DE  V  (  I  )  ,  PDEVI  I  ) 
,  UA.Si  ) 


7  0 

f  {1R*AT 

(  IH  1 

PL 

FORMAT 

(  [HO/ 

X.  w  *  ,  R  TmC  OMf'u  TE  L  Y 

VAL  UC  S 

AND  RESIDUALS 

(CONTINUED) 

9( 

con )  . 

j( 

'.0  TC 

l  fit; 

1  C( 

wR  1  1  F 

(  fc .  IK 

) 

!  1  ( 

FORMAT 

(  1 HC  / 

I  Owx  ,  U-PF  ot  E  NT  / 

lux,  SHI  ABEL-  1 S  X ,  1 HY 

,  1  S  «  ,  ,HT| 

L 

1  *>  *  * 

2H*2 ,  ! 

3  x  ,  7m  Y  L  *  L  .  t  i  1  X 

,  2 ( 6HY 

Ct V . ,  1  1  x  •  /  ) 

pH  I  3C 

I  «  1  , 

N 

Ml  « 

1  f' 

- 3  - 


F  12  -  FLOAT ( I ) / 5 - 0 

IF  (LIMES  .  L  T .  40)  GO  TO  115 

««ITE  (6,  20)  (TfTLF(K),  K  *  l,  16).  IPAGE 

IPAGE  *  I  PAGE  *  I 

wr  .  (6,  80) 

W  ITE  (6,  110) 

L  MfS  *  0 

1)5  WRITE  16,  120)  IV(I.K),  K  ■  1,  5),  YCll),  YOEVII),  POEV(I) 

120  FORMAT  I  IH  ,  9*,  t 4 ,  6(3X,  F14.S)  ) 

LINES  *  LINES  ♦  l 

IF  i  F I  1  „EQ.  F12)  WRITE  (6,  70) 

130  CONTINUE 
GO  TO  180 
1 A C  WRITE  (6,  150.' 

150  FORMAT  (IHO/  121X,  7MPERCENT  /  2 X ,  4HLABEL.  15X,  IH*,  15X,  2HXI, 

L  15X,  2 HX 2 *  l  5  X ,  2HX3,  13X,  Thy  CALC.,  11X,  6HY  lEV.,  IIX, 

C  6HY  DEV.  /  ) 

CO  IT?  I  »  l,  N 
FI  1  «  1/5 

F I  2  *  FLOAT! I ) / 6.0 

If  (LINES  .LT.  40)  GO  TO  155 

WRITE  (6,  20)  (TITLE(K).  X  •  l,  '  IMAGE 

IMAGE  -  IMAGE  ♦  1 

WRITE  (6,  80) 

WRITE  (6,  150} 

LINES  »  0 

155  WRITE  16,  160)  !VU,M,  *  •  1.  6),  YC  (1)  »  YDEV(l),  POEV(l) 

160  FORMAT  ( IH  ,  IX,  2A4,  7<3Xf  F14.5)  ) 

LIMES  «  LINES  ♦  l 

IF  ( F 1 1  .EG.  F 12 )  WRITE  (6,  70) 

1?0  CONTINUE 

180  IF  (NOTE  -EQ.  1)  WRITE  <6,  190) 

190  FORMAT  IlHO/,  10X,  109HTME  ABOVE  SOLUTION  FOR  THIS  CASE  DOES  NOT 
C  GIVE  AN  ABSOLUTE  MINIMIZATION  OF  THE  SUM  OF  SQUARES  OF  Y  RESIDUAL 
C  S.  ) 

RETURN 
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10.  ABSTRACT  Description  and  listing  of  an  air -FORTRAN  REY  WORDS 


IV  program  that  m^kes  least-squares  deter¬ 
minations  of  the  paranerers  of  any  of  five 
mathematical  functions  selected  by  the 
user,  given  a  set  of  observations  on  the 
'ependenr  and  independent  variables  of 
interest  (up  to  200  data  points  per  curve) . 
The  functions  available  are  those  most 
commonly  used  in  developing  cost  estimat¬ 
ing  relationships:  line,  parabola,  power, 
asymptotic-power,  and  exponential.  Up  to 
three  independent  variables  may  be  used 
for  the  line  and  power  functions.  The 
Y-intercept  may  be  specified  for  the 
line,  parabola,  or  asymptotic-power 
functions.  The  program  is  designed  to 
be  user-oriented  and  easily  workable 
rather  than  to  emphasize  computational 
efficiency  Exact  and  unique  solutions 
for  the  xine  and  parabolic  functions  are 
obtained  by  standard  algebraic  methods. 
Since  the  other  three  choices  are  not 
linear  in  all  parameters,  they  are  solved 
iteratively:  the  power  and  exponential  by 
a  modified  Gauss-Newton  iteration,  starting 
fron  the  exact  logarithmic  solution  (as 
described  in  RM~4y?9~PR),  and  the  asymp¬ 
totic-power  function  by  a  special  iterative 
method  described  in  die  RM.  In  addition  t: 
directions  for  pros>,ram  use,  the  RM  include, 
a  discussion  of  the  characteristics  of  the 
functions  and  mathematical  considerations 
involved  in  nonlinear  least-squares 
solutions . 
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